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Abstract 

Sets  of  yield  tables  that  show  probable  results  of  various  manage- 
ment alternatives  can  be  valuable  tools  for  decisionmaking,  es- 
pecially when  they  can  be  made  available  quickly  and  at  relatively 
low  cost.  Such  tables  can  be  obtained  with  data  from  temporary 
plots  and  the  computer  programs  presented. 

Key  words:  Stand  yield  tables,  timber  management,  simulation, 
managed  stands,  Pinus  ponder os a. 


USDA  Forest  Service 
Research  Paper  RM-79 


December  1971 


FIELD  AND  COMPUTER  PROCEDURES  FOR 
MANAGED-STAND  YIELD  TABLES 


Clifford  A.  Myers,  Principal  Mensurationist 
Rocky  Mountain  Forest  and  Range  Experiment  Station 


Central  headquarters  maintained  in  cooperation  with  Colorado  State 
University  at  Fort  Collins. 


Contents 

Page 

Introduction   1 

Uses  of  Yield  Tables   1 

General  Description  of  Methods   2 

Information  Needed   3 

1.  Stocking  After  Cutting   3 

2.  Description  of  Unthinned,  Young  Stands   5 

3.  Diameter  Increase  from  Growth   5 

4.  Diameter  Increase  from  Thinning   6 

5.  Height  of  Dominants  and  Codominants   6 

6.  Height  Increase  from  Thinning   6 

7.  Noncatastrophic  Mortality  

8.  Stand  Volume  Equation  

9.  Volume  Conversion  Factors  

Description  of  Program  PONYLD   8 

An  Application  of  PONYLD   1 1 

Modification  of  PONYLD   12 

Description  of  Program  PONCHK   I4 

Literature  Cited   I5 

Appendices: 

1.  Listing  of  Program  PONYLD   16 

2.  Output  of  PONYLD   19 

3.  Listing  and  Output  of  Program  PONCHK   22 


Field  and  Computer  Procedures  for 
Managed-Stand  Yield  Tables 


Clifford  A.  Myers 


Introduction 

Three  publications  have  presented  proce- 
dures for  computation  of  yield  tables  for  man- 
aged, even-aged  stands.  Two  of  these  described 
field  procedures  to  be  used  on  temporary  plots 
and  analysis  of  the  data  obtained  (Myers  1966, 
1967);  the  third  presented  a  computer  program 
that  calculated  and  printed  yield  tables,  using 
functions  described  in  the  other  papers  (Myers 
and  Godsey  1968). 

Field  procedures  and  computer  programs 
described  in  this  Paper  are  revisions  and  im- 
provements of  material  in  the  publications  cited 
above.  Some  modifications  are  in  response 
to  requests  by  users  for  applications  other  than 
those  intended  originally.  Other  changes  are 
the  result  of  increases  in  available  data  and 
greater  insight  obtained  by  applying  the  pro- 
cedures to  several  tree  species.  They  produce 
more  realistic  simulations  of  (1)  stand  growth, 

(2)  response  to  intermediate  cuttings,  and  (3) 
reproduction  cutting  by  any  of  the  even-aged 
systems.  Improved  procedures  are  used  for 
computations  such  as  periodic  mortality  and 
changes  in  average  height  and  diameter  with 
thinning. 

The  material  presented  includes  descriptions 
of:  (1)  field  measurements  needed  to  produce 
yield  tables  for  managed  stands,  (2)  relation- 
ships to  be  obtained  from  the  measurements, 

(3)  a  computer  program  written  in  FORTRAN 
IV  that  computes  and  prints  yield  tables,  (4) 
an  example  of  what  the  program  can  produce, 
and  (5)  a  computer  program  useful  in  generating 
certain  types  of  data.  Tabular  and  functional 
relationships  included  apply  to  ponderosa  pine 
(Pinus  ponderosa  Laws.)  in  the  Black  Hills  of 
South  Dakota  and  Wyoming.  Instructions  for 
adapting  the  procedures  and  programs  to  other 
species  or  forest  regions  are  given. 


Uses  of  Yield  Tables 

Yield  tables  for  managed  stands  are  essential 
guides  for  forest  managers.  They  report  probable 
wood  yields  that  result  from  specified  com- 
binations of  such  factors  as  site  quality,  utili- 
zation standards,  and  frequency  and  intensity 
of  thinning.  They  are,  therefore,  the  basis  for 
timber  management  planning.  They  also  pro- 
vide an  important  part  of  the  information  needed 
for  determining  the  influence  of  timber  treat- 
ments on  all  forest  resources. 

Yield  tables  for  a  species  are  useful  re- 
gardless of  the  current  level  of  management 
reached  in  forests  of  that  species.  Well-managed 
forests  can  benefit  from  refinements  in  opera- 
tions that  are  guided  by  comparisons  of  actual 
conditions  with  a  good  standard.  Where  con- 
version to  managed  stands  is  underway,  yield 
tables  provide  goals  toward  which  conversion 
can  be  directed.  Managers  can  better  examine 
alternatives  and  reach  decisions  concerning 
goals  through  computerized  simulation  of 
forestry  activities,  a  procedure  that  uses  yield 
table  computation  in  the  set  of  mathematical 
operations  (Myers  1968).  Several  factors  influ- 
ence the  ways  that  yield  tables  should  be 
produced  and  the  types  of  data  used  to  com- 
pute them.  These  factors  include: 

1.  Yield  tables  for  managed  stands  are  needed 
even  for  regions  and  species  where  managed 
stands  do  not  yet  exist.  The  term  "man- 
aged," as  used  here,  means  control  of  stand 
density  throughout  the  life  of  the  stand. 

2.  No  single  yield  table  can  be  adopted  as 
standard  for  a  species  or  region.  Applied 
forest  management  requires  the  use  of  various 
practices  that  can  and  will  differ  from  one 
forest  to  another,  and  must  appear  as  vari- 
ables in  yield  table  computation. 


-  1  - 


3.  A  manager  should  not  be  restricted  to  only 
one  yield  table  per  working  group,  or  series 
of  stands  managed  under  the  same  silvi- 
cultural  system.  He  must  have  the  oppor- 
tunity to  examine  the  probable  results  of 
his  operations,  to  make  necessary  changes 
in  the  management  of  any  of  his  resources, 
and  to  study  the  effects  of  these  changes 
before  money  is  spent  on  them. 

The  first  item  listed  above  appears  to  call 
for  the  use  of  temporary  plots  to  obtain  needed 
relationships.  The  next  two  sections  of  this 
Paper  describe  how  data  from  temporary  plots, 
permanent  plots,  or  combinations  of  both  can 
produce  the  relationships  needed  for  yield  table 
construction. 

The  second  and  third  items  above  relate  to 
a  disadvantage  of  normal  and  empirical  yield 
tables.  One  table  can  report  probable  yields 
for  only  one  combination  of  stand  characteristics 
and  management  objectives.  This  is  not  an 
effective  basis  for  comparing  alternatives  and 
making  decisions.  Since  a  group  of  yield  tables 
can  show  various  combinations  of  characteristics 
and  objectives,  a  means  of  producing  many 
yield  tables  at  low  cost  can  be  useful  to  a 
forest  manager.  Computer  programs  for  doing 
this  are  described  below  and  listed  in  the 
appendices.  Once  the  needed  relationships 
between  stand  variables  have  been  established, 
a  manager  can  examine  the  probable  outcomes 
of  many  possible  variations  in  management. 
There  is  no  need  to  delay  decisions  or  to  specu- 
late on  what  may  happen  if  a  condition  or 
procedure  changes.  Large  numbers  of  tables, 
each  based  on  a  specified  set  of  alternatives, 
can  be  computed  and  printed  at  a  cost  of  a 
few  cents  each. 


General  Description  of  Methods 

Nine  items  of  information,  described  in  the 
next  section,  are  needed  to  compute  yield  tables 
for  even-aged,  managed  stands.  The  first  item, 
the  residual  stand  desired  after  each  cutting,  is 
based  on  all  pertinent  information  available 
and  may  apply  to  several  species  and  regions. 
Data  for  the  second  item,  appearance  of  un- 
thinned  young  stands,  are  obtained  from  a  few 
measurements  made  on  temporary  plots.  The 
other  seven  items  are  based  on  data  from 
temporary  and /or  permanent  "growth-predic- 
tion plots"  that  are  measured  in  detail.  Analy- 
ses of  plot  data  described  below  depend  on  the 
availability  of  tree  volume  equations  and  site 
index  curves  applicable  to  the  species  and 


locality  (Bruce  and  Schumacher  1950,  Chapman 
and  Meyer  1949).  Procedures  for  estimation  of 
site  index  from  soil  and  topography  may  be 
used  where  stand  density  indicates  that  the 
conventional  height-age  relationship  will  give 
inaccurate  results  (Myers  and  Van  Deusen  1960). 

Growth-prediction  plots  are  temporary  or 
permanent  plots  placed  in  carefully  selected, 
even-aged  stands.  They  should  sample  numer- 
ous combinations  of  site  index  and  stand  density, 
including  extreme  values  of  possible  interest. 
The  sample  should  contain  a  full  representation 
of  ages  or  average  diameters,  at  least  to  the 
longest  rotations  for  which  yield  tables  may  be 
prepared.  Plots  must  also  conform  to  the  usual 
requirements  as  to  uniformity  of  site  quality, 
range  of  tree  sizes,  and  stand  density  across 
a  plot.  They  must  be  large  enough  to  provide 
an  adequate  representation  of  the  distribution 
of  trees  by  diameter  classes.  The  stands  must 
not  have  been  partially  cut  or  otherwise  dis- 
turbed during  the  period  for  which  growth  is 
measured.  They  must  not  be  diseased  or 
infested  by  insects  to  the  extent  that  growth 
is  affected,  unless  this  is  to  be  a  variable 
in  yield  table  construction. 

It  may  sometimes  be  desirable  or  necessary 
to  combine  data  from  temporary  and  permanent 
plots.  Or,  functions  can  be  derived  from  tem- 
porary plot  data  with  permanent  plots  used  as 
benchmarks  to  verify  accuracy  of  the  predictions. 
The  discussion  that  follows  assumes  that  tem- 
porary plots  will  provide  most  or  all  the  data 
needed. 

Most  forest  regions  contain  stands  of  various 
ages  and  with  densities  in  the  range  of  possible 
interest.  Such  stands  provide  a  readily  available 
store  of  growth  information.  Data  from  tem- 
porary plots  in  these  stands  are  entirely  adequate 
for  yield  table  construction  and,  except  for 
measures  of  periodic  mortality,  do  not  differ 
in  plot  values  or  variability  from  data  from 
equivalent  permanent  plots  (De'court  1965,  Myers 
1966,  Vuokila  1965). 

Measurements  made  on  each  growth-predic- 
tion plot  must  include  the  following: 

1.  Plot  area. 

2.  Heights  and  ages  of  trees  suitable  for  site 
index  determination.  Alternatively,  soil  and 
topographic  measurements  useful  in  esti- 
mating site  index  are  obtained. 

3.  D.b.h.  of  each  tree,  to  0.1  inch. 

4.  Total  height  of  each  tree  to  1.0  foot,  or  of 
each  of  a  sample  if  there  are  many  trees  on 
the  plot.  If  height  is  sampled,  sufficient 
measurements  must  be  taken  to  construct 
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a  height-diameter  curve  for  the  plot.  A  good 
sample  of  the  heights  of  dominant  and  co- 
dominant  trees  is  also  needed  so  the  average 
height  of  such  trees  can  be  computed. 

5.  Crown  class  of  each  tree. 

6.  Total  ages  of  a  sample  of  dominant  and  co- 
dominant  trees.  This  measurement  may  be 
omitted  if  even-aged  status  is  not  in  doubt 
or  if  equivalent  ages  are  obtained  for  site 
index  determination. 

7.  Radial  wood  growth  during  the  period  of 
interest,  to  0.05  inch,  from  increment  borings 
at  breast  height  along  the  best  estimates  of 
average  radius.  The  period  selected  will  be 
the  unit  prediction  period  of  the  computer 
program  that  computes  and  prints  the  yield 
tables. 

8.  D.b.h.  outside  bark  of  each  tree  that  died 
during  the  period  selected  for  item  7.  If 
necessary,  diameter  inside  bark  is  measured 
and  then  converted  to  diameter  outside  bark. 

9.  Several  cut  or  leave  codes  for  each  tree, 
based  on  trial  markings  of  each  plot  to 
simulate  several  intensities  of  thinning.  This 
information  is  used  to  determine  increases 
in  diameter  and  height  due  to  thinning,  as 
described  in  items  4  and  6  of  the  section 
headed  Information  Needed. 

Other  field  measurements  will  be  needed 
if  local  experience  indicates  that  additional 
independent  variables  are  significant  in  the 
regression  equations  described  below. 

Plot  and  tree  data  computed  initially  from 
the  field  measurements  are  as  follows: 

1.  Site  index. 

2.  Past  d.b.h.  of  each  tree  from  present  d.b.h., 
radial  wood  growth,  and  periodic  bark  growth 
(Myers  and  Van  Deusen  1958). 

3.  Height  of  each  tree  for  which  the  actual 
height  was  not  measured.  A  height-diameter 
curve  or  function  for  each  plot  provides 
any  missing  heights. 


Measured  and  computed  items  that  describe 
the  present  stand  are  used  to  compute  the 
following  values  for  each  plot.  Amounts  per 
acre  are  computed,  where  applicable. 

1.  Number  of  live  trees. 

2.  Basal  area. 

3.  Average  d.b.h.,  computed  as  the  tree  of 
average  basal  area. 

4.  Average  height  of  dominant  and  codominant 
trees. 

5.  Average  age  of  main-stand  trees. 

6.  Cubic  feet,  from  ground  to  tip,  of  all  trees. 


This  and  the  following  item  are  obtained  from 
tree  volume  equations. 
7.  Volumes  in  any  other  units  of  interest. 

Past  diameters  of  live  trees,  diameters  of 
the  tallied  dead  trees,  and  present  stand  age 
provide  the  following  items  that  describe  the 
stand  at  the  beginning  of  the  prediction  period. 

1.  Number  of  trees. 

2.  Basal  area. 

3.  Average  d.b.h.,  computed  as  the  tree  of 
average  basal  area. 

4.  Average  age  of  main-stand  trees,  from  present 
age  and  the  length  of  the  prediction  period. 

Once  the  plot  measurements  have  been  ob- 
tained and  summarized,  each  item  described 
in  the  next  section  is  computed  as  one  or  more 
relationships  to  be  placed  in  the  computer 
program.  These  functions  are  the  species- 
specific  elements  of  the  programs,  and  are 
replaced  if  the  program  is  to  apply  to  other 
species  or  utilization  standards. 


Information  Needed 

This  section  contains  instructions  for  ob- 
taining the  nine  items  that  appear  as  one  or 
more  FORTRAN  statements  each  in  PONYLD, 
the  computer  program  that  produces  yield  tables 
(appendix  1).  Most  of  the  items  are  obtained 
by  regression  analysis  of  plot  values  described 
in  the  previous  section.  Data  for  the  first  two 
items  described  come  from  other  sources,  in- 
cluding appropriate  references  in  the  literature. 
Additional  information  of  field  and  computation 
procedures  can  be  found  in  standard  mensura- 
tion texts  (Bruce  and  Schumacher  1950,  Chap- 
man and  Meyer  1949). 

Tabulations  include  only  enough  entries 
to  explain  the  nature  of  the  information  needed. 
They  do  not  indicate  sample  sizes  or  desirable 
ranges  of  data. 

Program  statements  relating  to  one  or  more 
of  the  nine  items  must  be  changed  to  adapt 
PONYLD  to  other  species  or  conditions.  Sug- 
gested program  modifications  are  given  in  the 
section  headed  Modification  of  PONYLD. 


1.  Stocking  After  Cutting. 

Stand  density  to  be  left  after  each  cutting 
is  expressed  as  a  relationship  between  basal 
area  and  average  stand  diameter  (table  1). 
Basal  area  is  an  easily  measured  variable  that 
is  highly  correlated  with  growth  in  size  and 
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Table  1. — Basal  areas  after  intermediate  cutting  in  relation  to  average  stand  diameter 

growing  stock  level  80 


Average 

Basal 

Average 

Basal 

Average 

Basal 

nVC l agC 

Basa.1 

stand  d.b.h. 

area 

stand  d.b.h. 

area 

stand  d.b.h. 

area 

<3  t  an  H   H  h 

B.  TQ  3. 

after  cutting 

per 

after  cutting 

per 

after  cutting 

3ft"PT  Piit"t"inp 

per 

(Inches) 

acre 

(Inches) 

acre 

(Inches) 

acre 

(Inches) 

acre 

Sq.  ft. 

Sq.  ft. 

Sq.  ft. 

Sq.  ft 

2.0 

12.1 

4.0 

35.2 

6.0 

56.6 

8.0 

72.5 

2.1 

13.2 

4.1 

36.4 

6.1 

57.6 

8.1 

73.1 

2.2 

14.4 

4.2 

37  .6 

6  .2 

58.5 

8 . 2 

7"?  7 

2.3 

15.5 

4.3 

38.7 

6.3 

59.4 

8.3 

74.3 

2.4 

16.7 

4.4 

39.9 

6.4 

60.3 

8.4 

74.8 

2.5 

17.9 

4.5 

41.0 

6.5 

61.2 

8.5 

75.3 

2.6 

19.0 

4.6 

42.2 

6.6 

62.1 

8.6 

75.8 

2.7 

20.2 

4.7 

43.4 

6  .7 

62  .9 

8 . 7 

76 . 3 

2.8 

21.3 

4.8 

44.5 

6.8 

63.8 

8.8 

76.7 

2.9 

22.5 

4.9 

45.7 

6.9 

64.6 

8.9 

77.1 

3.0 

23.7 

5.0 

46.8 

7.0 

65.4 

9.0 

77.5 

3.1 

24.8 

5.1 

47.8 

7.1 

66.2 

9.1 

77.9 

3.2 

26.0 

5.2 

48.8 

7.2 

67.0 

9.2 

78.2 

3.3 

27.1 

5.3 

49.8 

7.3 

67.7 

9.3 

78.5 

3.4 

28.3 

5.4 

50.8 

7.4 

68.5 

9.4 

78.8 

3.5 

29.5 

5.5 

51.8 

7.5 

69.2 

9.5 

79.1 

3.6 

30.6 

5.6 

52.8 

7.6 

69.9 

9.6 

79.3 

3.7 

31.8 

5.7 

53.8 

7.7 

70.6 

9.7 

79.5 

3.8 

32.9 

5.8 

54.7 

7.8 

71.2 

9.8 

79.7 

3.9 

34.1 

5.9 

55.7 

7.9 

71.9 

9.9 

79.8 

10.0+  80.0 


volume.  The  relationships  fit  well  into  com- 
puter algorithms  and  provide  thinning  guides 
for  application  of  treatments  to  actual  stands. 

Results  of  thinning  studies  and  data  from 
temporary  plots  are  used  to  construct  a  graph 
of  desired  basal  area  over  average  stand  diameter 
for  local  average  site  quality.  "Best"  stand 
density  for  each  average  diameter  sampled  can 
be  based  on  such  criteria  as  production  in  cubic 
feet  and  probable  length  of  saw  log  rotations. 
Selected  basal  areas  are  plotted  over  corre- 
sponding stand  diameters,  and  a  curve  is  drawn 
through  the  points.  If  the  relationship  appears 
reasonable,  it  is  then  expressed  as  a  mathe- 
matical function. 

In  table  1,  basal  area  increases  with  diameter 
until  10.0  inches  diameter  is  reached,  and  re- 
mains constant  thereafter.  The  designation 
"growing  stock  level  80"  indicates  that  basal 
area  is  80.0  square  feet  when  diameter  is  10.0 
inches  or  larger,  regardless  of  what  basal  area 
may  be  at  lower  average  diameters. 


Desired  stand  density  will  vary  with  the 
objectives  of  management,  and  a  family  of 
basal  area-diameter  relationships  is  needed  (fig. 
1.).  The  original  single  curve  or  function  of 
basal  area  on  diameter  is  treated  as  a  guide 
curve  from  which  other  curves  can  be  pro- 
duced. Basal  areas  for  any  growing  stock 
level  can  be  computed  by  multiplying  the  guid- 
ing values  for  level  80  in  table  1  by  the  ratio 
level/80.  For  example,  basal  areas  for  level 
100  are  each  100/80  times  the  corresponding 
basal  areas  of  table  1.  The  level  designations 
are  the  variables  THIN,  DLEV,  and  DSTY  in 
PONYLD  (appendix  1). 

The  curves  of  figure  1  define  growing  stock 
goals  for  many  possible  management  objectives, 
including  the  possibility  of  holding  stand  density 
above  or  below  levels  that  give  the  best  saw 
log  or  roundwood  production.  Any  desired 
form  of  the  guide  curve  may  be  used  if  the 
appropriate  statements  of  PONCUT  are  modi- 
fied properly. 
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Relationships  shown  in  table  1  appear  as 
functions  for  level  80  in  program  PONYLD. 
Basal  areas  computed  from  these  functions 
are  multiplied  by  terms  that  include  the  desired 
growing  stock  level  (THIN  or  DLEV)  to  obtain 
values  for  other  growing  stock  levels.  Variables 
for  which  FORTRAN  statements  are  needed 
and  their  use,  are: 

a.  DBHP  -  to  find  a  d.b.h.  less  than  10.0  inches 
when  basal  area  is  known.  Three  equations 
for  DBHP  are  used  to  simplify  representation 
of  the  nonlinear  relationship  between  d.b.h. 
and  basal  area. 

b.  BREAK  and  BUST  -  to  compute  values  of 
basal  area  that  are  the  upper  limits  of  appli- 
cability of  the  first  two  equations  for  DBHP. 

c.  SQFT  -  to  find  basal  area  when  d.b.h.  is 
known.  Two  equations  represent  the  non- 
linear relationship  for  d.b.h.  less  than  10.0 
inches. 


2.  Description  of  Unthinned,  Young  Stands. 

Values  in  the  first  line  of  each  yield  table 
describe  stand  conditions  just  prior  to  initial 
thinning.  They  are  entered  directly  from  data 
cards  or  are  computed  from  the  data.  Users 
of  PONYLD  must,  therefore,  be  able  to  describe 
the  stands  that  do  or  should  exist  at  time  of 


initial  thinning.  In  a  yield  table  for  managed 
stands,  the  stand  density  and  related  average 
diameter  given  in  the  first  line  result  when 
stand  regeneration  and  subsequent  growth  and 
mortality  progress  as  planned. 

Numerous  unthinned,  young  stands  can  be 
examined  to  determine:  (1)  site  index,  (2) 
average  stand  d.b.h.,  (3)  number  of  trees  per 
acre,  and  (4)  stand  age.  In  selecting  plots, 
preference  should  be  given  to  those  stands 
that  represent  possible  regeneration  goals  for 
various  objectives  of  management.  For  each 
site  index  class,  average  d.b.h.  is  determined 
for  various  combinations  of  stand  age  and  num- 
ber of  trees  per  acre.  Any  influence  of  the 
overstory  must  be  included  where  shelterwood 
or  seed  tree  systems  are  used. 


3.  Diameter  Increase  From  Growth. 

Regression  analysis  of  data  obtained  on  the 
growth  prediction  plots  described  in  the  previous 
section  provides  an  equation  for  prediction  of 
future  average  stand  d.b.h.  In  terms  of  the 
plot  data  described  above,  present  average  d.b.h. 
of  ponderosa  pine  is  estimated  from  past  d.b.h., 
site  index,  and  past  basal  area  per  acre.  Addi- 
tional independent  variables  may  be  useful  for 
other  species  or  localities.  The  prediction 
period  is  determined  by  the  number  of  rings 
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measured  on  increment  cores.  The  equation 
for  ponderosa  pine  and  for  a  10-year  prediction 
period  appears  in  PONYLD  as  the  FORTRAN 
statement  for  DBHO. 

4.  Diameter  Increase  From  Thinning. 

Change  in  average  stand  diameter  caused 
by  intermediate  cuttings  can  be  estimated  from 
data  obtained  during  repeated  trial  marking  of 
the  plots  described  previously.  Numerous  inten- 
sities of  cutting  are  simulated  on  each  plot, 
and  each  tree  is  coded  as  cut  or  leave  at  each 
trial.  Since  the  dimensions  of  each  plot  tree 
are  known,  volumes  and  other  values  can  be 
determined  for  each  potential  residual  stand. 
Values  useful  in  computation  of  prediction 
equations  for  diameter  change  and  other  items 
include:  average  d.b.h.,  basal  area  per  acre, 
percentage  of  trees  retained,  average  height, 
volumes  per  acre  of  entire  stems  in  cubic  feet 
(total  cubic  feet),  and  volumes  in  other  units 
of  interest. 

In  PONYLD,  diameter  after  thinning  is 
estimated  from  diameter  before  thinning  and  the 
percentage  of  trees  to  be  retained.  Regression 
analysis  of  data  from  the  simulated  thinnings 
provide  functions  for  DBHE  andPDBHE.  These 
two  variables  represent  the  same  item,  diameter 
after  thinning.  DBHE  is  computed  directly 
if  the  estimated  percentage  of  trees  to  be  re- 
tained is  at  least  50  percent.  With  fewer  trees 
retained,  the  relationship  is  highly  nonlinear, 
so  PDBHE  is  computed  and  its  antilogarithm 
becomes  DBHE. 

Because  of  the  numerous  possible  com- 
binations of  initial  stand  density,  stocking  level, 
and  initial  stand  diameter,  trial  marking  of 
sample  plots  may  not  provide  all  the  information 
needed.  A  supplementary  procedure  is  described 
in  the  section  headed  Description  of  Program 
PONCHK  and  in  appendix  3. 


5.  Height  of  Dominants  and  Codominants. 

Average  heights  of  dominant  and  codomi- 
nant  trees  are  computed  from  data  obtained  on 
growth-prediction  plots  where  height  growth 
apparently  has  never  been  reduced  by  high  stand 
density.  Regression  analysis  of  the  data  provides 
functions  for  estimating  height  for  various 
combinations  of  site  index  and  stand  age  (table 
2).  These  functions  appear  as  statements  for 
HTSO  in  PONYLD.  Two  functions  are  used 
to  cover  the  range  of  possible  ages.  With  the 
procedure  used  to  account  for  changes  in  height 
from  thinning,  average  d.b.h.  cannot  be  used 
instead  of  site  index  and  age  to  predict  height. 


Table  2. — Average  height  (feet)  of  dominant 
and  codominant  tre*«  at  various 
ages,  Black  Hills  ponderosa  pine 
(as  computed  by  statements  for 
HTSO) 


Main  Site  index  class 

stand  age  

(Years)  40  50  60  70 


20 

8 

10 

1 1 
±  j 

1  A 
-L  D 

30 

12 

16 

21 

25 

40 

17 

22 

28 

34 

50 

21 

28 

36 

43 

60 

26 

33 

41 

49 

70 

30 

39 

47 

56 

80 

34 

43 

52 

62 

90 

37 

47 

57 

66 

100 

40 

50 

60 

71 

110 

43 

53 

64 

74 

120 

45 

56 

67 

77 

130 

47 

58 

69 

80 

140 

48 

60 

71 

82 

150 

50 

62 

73 

85 

160 

51 

63 

75 

87 

Table  2  is  similar  to  table  1  of  Technical 
Bulletin  630  (Meyer  1938).  This  suggests  that 
heights  from  good  site  index  curves  or  tables 
may  supplement  or  substitute  for  local  data, 
until  adequate  sampling  can  be  completed.  If 
site  index  data  are  used,  they  must  be  based 
on  the  same  crown  classes  as  the  equations 
for  stand  volume  in  cubic  feet,  described  below. 


6.  Height  Increase  From  Thinning. 

Changes  in  the  average  height  of  dominant 
and  codominant  trees  due  to  partial  cutting  are 
estimated  the  same  way  as  changes  in  average 
diameter.  Results  of  repeated  trial  markings 
provide  the  data  needed  to  compute  the  pre- 
diction function.  In  PONYLD,  the  variable 
ADDHT  is  the  computed  amount  of  change. 
The  percentage  of  trees  retained  is  the  only 
independent  variable  used.  At  each  cutting, 
the  current  value  of  ADDHT  is  added  to  height 
before  thinning,  HTSO,  to  obtain  height  after 
thinning,  HTST.  It  is  also  added  to  a  cumu- 
lative sum  of  changes,  HTCUM,  so  computed 
heights  before  thinning  will  show  the  effects 
of  past  treatment  as  well  as  of  age.  Computed 
values  of  ADDHT  are  small  because  average 
heights  and  their  changes  refer  to  dominant 
and  codominant  trees,  only. 
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As  with  diameter,  it  may  be  difficult  to 
sample  many  combinations  of  stocking  level, 
initial  stand  density,  and  initial  average  height. 
The  supplementary  procedure  described  in  the 
section  Description  of  Program  PONCHK  may 
be  useful  in  such  cases. 


7.  Noncatastrophic  Mortality. 

Reduction  in  numbers  of  trees  may  be 
important  in  unthinned  stands,  but  minor  and 
erratic  in  thinned  stands.  Such  is  the  case  with 
ponderosa  pine  in  the  Black  Hills.  A  prediction 
equation  for  mortality  could  not  be  computed 
for  thinned  stands  with  an  average  d.b.h.  of 
10.0  inches  or  larger.  Small  reductions  are 
made  in  PONYLD  for  thinned  and  unthinned 
stands  of  smaller  diameter. 

The  prediction  equation  in  PONYLD  for 
percentage  of  mortality  expressed  as  a  decimal 
(DIED)  contains  average  stand  d.b.h.  and  basal 
area  per  acre,  both  at  the  beginning  of  the 
period,  as  independent  variables.  Data  to  com- 
pute such  equations  can  be  obtained  from 
two  sources: 

a.  Permanent  plots  that  have  been  measured 
at  least  as  frequently  as  the  prediction  period 
to  be  used. 

b.  Growth-prediction  and  other  temporary  plots 
that  have  not  been  partially  cut  for  at  least 
a  number  of  years  equal  to  the  prediction 
period.  It  must  be  possible  to  estimate  the 
number  of  years  since  death  for  each  tree 
that  died  during  the  prediction  period. 


8.  Stand  Volume  Equation. 

Plot  tallies  of  tree  diameters  and  heights 
are  converted  to  volumes  per  acre  and  other 
stand  measures,  as  already  described.  Data 
from  growth  prediction  plots,  results  of  simu- 
lated thinnings,  and  other  available  stand  tallies 
are  used.  Volumes  are  computed  in  total  cubic 
feet  and  in  other  units  of  interest  with  appro- 
priate tree  volume  equations  (Myers  1964). 
Only  total  cubic  feet  will  be  used  to  compute 
stand  volume  equations.  Total  cubic  volume 
per  acre  from  ground  line  to  tip  of  all  trees 
more  than  4.5  feet  tall  is  the  only  volume  com- 
puted directly  by  PONYLD.  It  is,  therefore, 
the  only  unit  for  which  a  stand  volume  equa- 
tion will  be  needed. 

Two  forms  of  the  stand  volume  equation 
have  proven  useful.  They  are: 


V  =  a  +  b,  BH  +  b2  D  (1) 

V  =  (a  +  b,D2H  +  b2B)xN  (2) 

where 

V  =  gross  total  cubic  volume  per  acre. 
B  =  basal  area  per  acre  in  square  feet. 
H  =  average  height  of  dominant  and  codominant 
trees. 

D  =  average  stand  d.b.h.  in  inches. 
N  =  number  of  trees  per  acre. 

Stand  volume  equations  appear  as  state- 
ments to  compute  CUFT.  Two  statements  are 
used  because  the  relationship  is  not  linear 
over  the  ranges  of  BH  or  D  2  H  that  can  appear 
in  the  yield  tables. 

Stand  volume  equations  must  apply  over  a 
wide  range  of  stand  densities  and  average 
diameters.  Procedures  described  in  the  section 
Description  of  Program  PONCHK  may  be  useful 
in  providing  the  variety  of  conditions  needed 
for  regression  analysis  or  test  of  the  functions. 


9.  Volume  Conversion  Factors. 

Program  PONYLD  computes  volumes  in 
units  other  than  total  cubic  feet  of  bole  wood 
from  cubic  volumes  per  acre  and  appropriate 
conversion  factors.  Plot  data  that  produce 
stand  volume  equations  must,  therefore,  also 
provide  conversion  factors.  Data  from  growth 
prediction  plots,  trial  thinnings,  and  other 
sources  are  used  to  compute  plot  volumes  in 
various  units,  as  described  previously.  Mer- 
chantable cubic  feet,  board  feet,  square  feet 
of  veneer,  weight  per  acre,  or  other  units 
may  be  computed.  Then,  the  quantity  of  each 
unit  per  total  cubic  foot  is  determined  separately 
for  each  plot.  Selection  of  appropriate  units 
includes  choice  of  minimum  merchantable  top 
diameters  and  diameters  at  breast  height.  The 
second  column,  below,  shows  some  of  the  ratios 
used  to  obtain  equations  for  FCTR.  The  third 
column  shows  ratios  used  for  PROD. 


Average  Merchantable  cubic      Board  feet 

stand  d.b.h.  feet  t  total  t  total 

(Inches)  cubic  feet  cubic  feet 


5 

1 

0.355 

6 

0 

.552 

6 

9 

.725 

8 

3 

.860 

0 

99 

9 

1 

.901 

1 

55 

10 

3 

.931 

2 

38 

19 

0 

.962 

5 

33 

23 

4 

.969 

5 

88 
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Ratios  need  not  be  computed  for  the  smallest 
average  stand  diameters  because  the  factors 
vary  greatly  in  such  stands.  This  has  no 
important  effect  on  yield  table  construction. 
Merchantable  material  will  not  be  a  part  of 
thinnings  until  stand  diameter  is  large  enough 
for  conversion  factors  to  be  reliable. 

Regression  analysis  provides  functions  to 
compute  conversion  factors  from  other  stand 
variables.  FCTR  can  be  estimated  from  average 
stand  d.b.h.  Estimates  for  PROD  for  ponderosa 
pine  are  improved  if  basal  area  per  acre  is 
also  included  in  the  equation.  More  than 
one  equation  for  FCTR  and  PROD  appear  in 
PONYLD  so  the  relationships  can  be  expressed 
by  simple  linear  functions  over  a  wide  range 
of  d.b.h. 

As  with  stand  volume  equations,  conver- 
sion factors  must  apply  over  a  wide  range  of 
stand  conditions.  The  description  of  program 
PONCHK,  below,  presents  a  solution  to  this 
problem. 


Description  of  Program  PONYLD 

Program  PONYLD  consists  of  a  main  pro- 
gram and  two  subroutine  subprograms  written 
in  standard  FORTRAN  IV.  The  main  program 


reads  the  data  cards,  performs  most  compu- 
tations, and  writes  the  yield  tables.  Subroutine 
PONCUT  determines  the  new  average  stand 
diameter  after  cutting  to  the  specified  growing 
stock  level.  Subroutine  PONVOL  computes 
volumes  in  total  cubic  feet  per  acre,  and  factors 
to  convert  these  volumes  to  other  units.  Opera- 
tions performed  by  each  routine  are  indicated 
by  comment  statements  in  the  source  program 
(appendix  1). 

Initial  operations  are  the  reading  of  the  five 
data  cards  described  in  the  tabulation  of  the 
order  and  contents  of  the  data  deck.  The  first 
two  cards  enter  values  that  will  not  change 
during  a  computer  run.  One  variable,  NTSTS, 
controls  the  number  of  sets  of  yield  tables 
to  be  produced,  one  for  each  group  of  cards 
of  types  3,  4,  and  5.  Other  variables  determine 
which  volumes  from  intermediate  cuts  will 
be  included  in  the  total  yields  printed  at  the 
bottom  of  each  table,  and  the  basis  of  the 
growing  stock  levels.  Three  data  cards  are 
read  for  each  test.  They  enter  the  number  of 
tables  to  be  produced  in  a  set,  initial  stand  con- 
ditions, and  controls  on  the  operations  to  be 
performed.  Entry  of  a  negative  or  zero  value 
for  any  variable  on  card  types  1  to  4  or  for 
REGN(l)  on  card  type  5  will  cause  the  printing 
of  an  error  message  and  termination  of  the  job. 


Order  and  Contents  of  the  Data  Deck  for  Program  PONYLD 


Card      Number  of  Variable 

type  cards         name         Columns       Format       Description  of  variable 


1        NTSTS  1-4 


GIDE  5-8 


1        COMCU  1-8 


COMBF  9-16 


1  per  JCYCL  1-4 
test 


14 


F4.0 


F8.3 


F8.3 


14 


Number  of  tests  per 
batch.  The  number  of 
sets  of  yield  tables  to  be 
produced. 

Base  level  of  set  of  grow- 
ing stock  levels,  as  the 
80.  in  table  1  and  figure  1. 
Minimum  cut  in  mer- 
chantable cubic  feet  to  be 
included  in  total  yields. 
Must  be  at  least  1.0. 
Minimum  cut  in  board 
feet  to  be  included  in 
total  yields.  Must  be  at 
least  1.0. 

Interval  between  inter- 
mediate cuts.  A  multiple 
of  RINT. 
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Card      Number  of  Variable 

type  cards        name         Columns       Format       Description  of  variable 


MIX 


1  per  AGEO 
test 


DBHO 


DENO 


DSTY 


RINT 


SITE 
THIN 


lper  REGN(l) 
test 


5-8 


1-8 


9-16 


17-24 


25-32 


33-40 


41-48 
49-56 


1-8 


VLLV(l)  9-16 


INVL(l)  17-24 


14 


F8.3 


F8.3 


F8.3 


F8.3 


F8.3 


F8.3 
F8.3 


F8.3 


F8.3 


F8.3 


Number  of  stocking  lev- 
els or  values  of  DLEV 
to  be  examined  in  one 
test. 

Initial  age  to  be  shown 
in  a  yield  table.  Stand 
age  when  first  thinning 
occurs. 

Average  stand  d.b.h.  just 
prior  to  initial  thinning 
at  stand  age  AGEO. 
Number  of  trees  per 
acre  just  prior  to  initial 
thinning  at  stand  age 
AGEO. 

Lowest  growing  stock 
level  for  intermediate 
cuts  after  initial  thinning. 
Level  will  increase  by  10 
as  many  times  as  speci- 
fied by  MIX  on  card  type 
3. 

Number  of  years  for 
which  a  growth  equation 
makes  one  projection. 
Value  is  10.0  for  state- 
ments in  appendix  1. 
Site  index  for  the  species. 
Growing  stock  level  for 
initial  thinning  at  age 
AGEO.  May  equal  DLEV. 
Stand  age  at  which  first 
regeneration  cut  will  oc- 
cur. Must  never  be  zero 
or  blank,  as  REGN  (1)  is 
rotation  length  for  clear- 
cutting. 

Percentage  of  previous 
DLEV  to  be  left  at  age 
REGN(l).  Will  be  zero 
with  clearcutting. 
New  interval  between 
cuts  in  effect  after  age 
REGN(l).  Will  be  zero 
with  clearcutting. 
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REGN(2)  25-32 


F8.3 


VLLV(2)  33-40 


INVL(2) 
REGN(3) 


41-48 
49-56 


F8.3 


F8.3 


F8.3 


Stand  age  at  which  second 
regeneration  cut,  if  any, 
will  occur.  Removal  of 
seed  trees  or  second  cut 
of  shelterwood. 
Percentage  of  previous 
DLEV  (including  effect 
of  VLLV(l))  to  be  left 
at  age  REGN(2).  May 
be  zero. 

New  interval  between 
cuts  in  effect  after  age 
REGN(2).  May  be  zero. 
Stand  age  at  which  third 
regeneration  cut,  if  any, 
will  occur.  Final  cut  of 
3-cut  shelterwood. 


All  operations  from  statement  30  to  state- 
ment 160  are  performed  for  each  table.  The 
tables  of  a  set  differ  by  the  growing  stock  level 
to  be  left  after  the  second  or  later  thinnings 
(DLEV).  The  first  table  of  each  set  is  com- 
puted with  the  subsequent  level  entered  as 
DSTY  on  card  type  4.  Each  following  table 
will  be  computed  with  second  and  later  thinnings 
cut  to  a  level  10  higher  than  that  of  the  pre- 
vious table.  Operations  performed  for  each 
table  are: 

1.  Computation  of  height  and  volumes  just 
prior  to  initial  thinning. 

2.  Partial  cutting  to  the  growing  stock  level 
specified  by  THIN  for  initial  thinning  or  by 
DLEV  for  subsequent  cutting.  Cutting  will 
not  be  simulated  if  the  stand  is  already 
below  the  basal  area  specified  by  THIN  or 
DLEV. 

3.  Computation  of  post-cutting  volumes  and 
other  values. 

4.  Printing  of  table  headings  the  first  time 
through  the  loop  and,  each  time  through 
the  loop,  printing  of  values  appropriate  for 
the  stand  age. 

5.  Projection  of  diameter,  height,  and  stand 
density  for  one  or  more  periods  until  the 
next  intermediate  cut  is  scheduled.  Stand 
volumes  and  other  values  are  computed  and 
printed  at  ages  when  no  cutting  is  scheduled. 

6.  Repetition  of  steps  2  to  5  until  stand  age  at 
time  of  initial  regeneration  cutting  is  reached. 

7.  Redefinition  of  the  growing  stock  to  be  left 
after  cutting  and  the  interval  between  re- 
generation cuttings.  Operations  are  guided 
by  entries  on  card  type  5.  Clearcutting 


calls  for  only  one  entry,  stand  age  at  time 
of  removal  or  REGN(l).  Seed  tree  cutting 
requires  entry  of  values  for  all  items  on 
card  type  5  to  and  including  REGN(2)  which 
is  the  age  when  seed  trees  are  to  be  re- 
moved. The  interval  between  regeneration 
cutting,  REGN(l),  and  removal  of  seed  trees 
is  INVL(l).  Shelterwood  cuttings  are  con- 
trolled similarly  except  that  up  to  three 
cuts  are  possible.  Stand  age  at  final  cut 
will  be  REGN(2)  for  two-cut  shelterwood 
and  REGN(3)  for  three-cut  shelterwood. 

8.  Regeneration  cuts  are  accomplished  by  repeti- 
tion of  steps  2  to  7,  until  the  age  of  final 
cut  is  reached. 

9.  Printing  of  column  totals  for  volumes  re- 
moved. Volumes  less  than  COMCU  and 
COMBF  on  card  type  2  will  not  be  included 
in  the  totals  so  that  commercial  yields  may 
be  compared,  if  desired.  Actual  column 
totals  may  be  obtained  by  entering  values 
of  1.0  for  COMCU  and  COMBF  on  card 
type  2. 

Subroutine  PONCUT  computes  average 
stand  diameter  after  cutting  from  diameter  before 
cutting  and  the  percentage  of  trees  retained. 
The  percentage  of  trees  retained  is  needed 
as  an  independent  variable,  but  is  itself  an 
unknown.  Successive  percentages  of  trees  are 
therefore  tested  until  d.b.h.  after  thinning, 
number  of  trees  retained,  and  residual  basal 
area  agree  with  the  diameter  and  basal  area 
combination  called  for  by  THIN  or  DLEV. 

These  combinations,  shown  in  table  1,  ap- 
pear in  PONCUT  as  statements  for  DBHP  in 
the  second  loop  and  as  REST  in  the  first  loop. 
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Two  loops  provide  for  the  increasing  and  con- 
stant segments  of  each  growing  stock  level 
line  (fig.  1).  Limiting  d.b.h.  for  selection 
of  loops  is  10.0  inches  minus  the  average 
change  expected  over  a  range  of  cuttings  in 
stands  just  under  10.0  inches  d.b.h. 

Subroutine  PONVOL  computes  total  cubic 
feet  per  acre  and  factors  to  convert  this  to 
other  units.  Conversions  to  merchantable  cubic 
feet  and  to  board  feet  are  shown  in  the  listing 
in  appendix  1.  Utilization  standards  for  these 
units  are  given  in  the  comment  statements 
of  PONVOL.  Conversions  to  other  units  or 
utilization  standards  may  supplement  or  replace 
those  already  in  PONVOL.  Additional  units 
might  be  square  feet  of  veneer  or  wood  weight 
in  pound?  (Myers  1960). 

Program  PONYLD  should  run  with  little 
or  no  modification  on  any  computer  that  ac- 
cepts FORTRAN  IV,  has  a  minimum  of  32K 
words  of  storage,  and  has  two  input/output 
devices  (unit  5  for  program  and  data  deck 
input  and  unit  6  for  printed  output).  Changes 
to  adapt  the  program  to  other  tree  species 
or  utilization  standards  and  for  additional  com- 
putations are  described  in  the  section  headed 
Modification  of  PONYLD. 


An  Application  of  PONYLD 

The  problem  described  below  demonstrates 
the  computations  made  by  PONYLD  and  the 
printed  results  obtained.  It  illustrates  some 
of  the  questions  that  may  be  asked  and  the 
information  that  will  be  provided.  The  example 
also  serves  as  a  test  problem  for  use  in  adapt- 
ing the  source  program  to  locally  available 
computing  facilities. 

A  forest  manager  wishes  to  determine  the 
intensity  of  thinning  that  will  maximize  volume 
production  in  board  feet  in  stands  of  site  index 
70.  Length  of  the  cutting  cycle  has  not  been 
standardized,  but  will  be  20  years  for  this  test. 
He  also  wants  to  compare  yields  from  two-cut 
and  three-cut  shelterwood,  both  with  the  final 
removal  cut  scheduled  for  stand  age  130  years 
and  considering  the  current  crop,  only.  Alter- 
natives calling  for  more  than  one  precommercial 
thinning  are  unacceptable.  Minimum  com- 
mercial volumes  per  acre  are  320  cubic  feet 
to  a  4-inch  top  and  1,500  board  feet.  The 
manager  expects  that  his  procedure  for  re- 
generation cuts  will  result  in  a  new  stand 
that  contains  950  trees  per  acre  by  age  30, 
with  an  average  diameter  of  4.8  inches. 

The  data  deck  consists  of  the  32  cards 
shown  in  figure  2.    The  first  card  calls  for 
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10  tests,  so  there  will  be  10  sets  of  one  each 
of  card  types  3,  4,  and  5.  The  second  card 
enters  minimum  merchantable  volume  limits. 
The  third  card,  a  type  3  card  that  is  repeated 
nine  more  times,  enters  the  20-year  cutting 
cycle  and  information  that  five  subsequent 
thinning  levels  are  to  be  examined  in  each 
test.  With  10  tests  of  five  thinning  levels 
each,  the  computer  will  produce  50  yield  tables. 
For  brevity,  only  a  few  of  them  are  repro- 
duced in  appendix  2. 

The  fourth  card  of  the  data  deck,  a  type 
4  card  that  follows  the  type  3  card,  differs 
from  other  type  4  cards  only  in  the  thinning 
level  that  controls  initial  thinnings.  For  each 
of  the  two  shelterwood  cuttings,  levels  80, 
90,  100,  110,  and  120  will  each  be  imposed  at 
age  30. 

The  fifth  card  is  the  type  5  card  for  two- 
cut  shelterwood.  A  removal  cut  that  retains 
half  the  basal  area  specified  for  the  subsequent 
thinning  level  is  scheduled  for  age  110.  The 
final  cut  will  be  made  at  stand  age  130.  The 
type  5  card  is  the  same  for  the  first  five  tests 
to  be  performed.  Card  number  20  of  the  data 
deck  is  also  a  type  5  card.  This  card  and  the 
other  four  type  5  cards  specify  three-cut  shelter- 
wood  with  cuts  at  ages  90,  110,  and  130.  The 
first  removal  cut  leaves  two-thirds  of  the  basal 
area  of  the  subsequent  thinning  level.  The 
second  removal  cut  leaves  half  the  basal  area 
left  after  the  first  removal  cut. 

Yield  tables  produced  by  PONYLD,  a  few 
of  which  are  reproduced  in  appendix  2,  can 
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assist  in  decisionmaking  in  many  ways.  Money 
yields  and  rates  earned  can  be  computed  by 
applying  thinning  costs  and  stumpage  values 
to  the  volumes  given  in  the  tables.  Stand 
ages  at  culmination  of  mean  annual  increment, 
and  rates  earned  assist  in  the  selection  of 
rotations. 

For  the  situation  described  above,  yields 
and  numbers  of  precommercial  thinnings  are  of 
greatest  immediate  interest.  These  items  are 
summarized  in  tables  3  and  4  for  the  50  yield 
tables  produced.  Combinations  of  low  initial 
and  low  subsequent  growing  stock  levels  or  of 
high  initial  and  intermediate  subsequent  levels 
produce  the  greatest  volumes  with  one  pre- 
commercial thinning.  Additional  comparisons 
should  be  made  to  include  such  factors  as 
probaDle  thinning  costs,  cubic  yields  from  thin- 
nings not  commercial  for  board  feet,  and  the 
average  size  of  tree  produced.  As  expected, 
the  current  crop  produces  more  board  feet  in 
130  years  if  cut  by  two-cut  shelterwood  than  if 
by  three-cut  shelterwood.  The  latter  treatment 
may,  however,  get  the  next  crop  off  to  an 
earlier  start. 

The  manager  may  now  wish  to  vary  the 
interval  between  cuts  and  other  variables  subject 
to  his  control.  The  production  of  many  more 
yield  tables  is  entirely  practical.  The  50  tables 
already  examined  required  only  17.9  seconds  of 
central  processor  time  on  a  CDC  6400  computer. 


Table  4. — Number  of  precommercial  thinnings  if 
each  of  the  50  combinations  of 
initial  and  subsequent  thinning 
levels  is  executed  as  specified  by 
the  data  deck  (both  types  of  cutting 
gave  the  same  results) 
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Modification  of  PONYLD 

PONYLD  can  be  modified  in  many  ways  to 
answer  questions  of  the  "what  would  happen 
if  I  did  this"  type.  Several  possible  modifica- 
tions are  listed  here. 


To  Change  Species: 


Table  3. — Yields  in  board  feet,  including 
commercial  thinnings,  of  the 
50  combinations  of  initial  and 
subsequent  thinning  levels. 
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PONYLD  can  be  converted  to  apply  to  a 
different  species  than  the  ponderosa  pine  used 
in  the  example.  This  change  requires  new 
statements  to  compute  HTSO,  ADDHT,  DBHO, 
and  DENO  in  the  main  program.  Replacement 
of  statements  for  CUFT,  FCTR,  and  PROD  are 
needed  in  PONVOL  and  for  DBHE  and  PDBHE 
in  PONCUT.  Another  definition  of  growing 
stock  than  that  shown  in  table  1  and  figure  1 
requires  new  statements  for  BREAK,  BUST, 
DBHP,  and  SQFT.  A  new  base  value  (GIDE) 
can  be  entered  on  card  type  1  if  a  level  other 
than  80  is  the  base  of  the  new  growing  stock 
levels.  Ways  of  obtaining  function  for  these 
variables  are  given  in  the  section  headed  Infor- 
mation Needed. 


To  Include  Losses: 

Modifications  can  be  made  to  include  the 
effects  of  insects  or  diseases  on  tree  growth 
and  stand  density.  How  this  may  be  done  for 
dwarf  mistletoe  is  explained  elsewhere  (Myers 
etal.  1971). 
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To  Add  Other  Measures: 

Quantities  other  than  potential  wood  yields 
can  be  computed.  For  example,  stand  measure- 
ments used  to  compute  volumes  can  also  be 
used  to  compute  the  amount  of  slash  that 
would  be  produced  by  cutting  the  quantities 
specified  in  the  tables. 


To  Study  Actual  Stands: 

Analysis  of  actual,  thinned  stands  has  been 
a  frequently  used  modification  of  the  original 
version  of  PONYLD.  In  this  case,  stand  condi- 
tions after  initial  cutting  are  known,  and  informa- 
tion on  possible  response  to  this  and  future 
treatments  is  desired.  Alternatively,  the  un- 
marked component  of  the  stand  is  known  and 
the  purpose  of  analysis  is  to  determine  if  treat- 
ment is  justified.  Modification  may  be  accom- 
plished with  the  following  changes  in  the  main 
program: 

1.  Read  AGEO,  DBHO,  and  DENO  as  usual, 
to  describe  the  actual  stand  just  prior  to 
thinning. 

2.  Add  a  READ  statement  after  the  READ 
statement  for  card  type  4  to  enter  actual 
values  of  DBHT,  DENT,  HTST,  and  BAST. 
These  describe  the  stand  immediately  after 
initial  thinning. 

3.  Change  the  DIMENSION  statement  from 
VAR(IO)  to  VAR(14). 

4.  Follow  the  READ  card  of  step  2  with  five 
statements: 

VAR(ll)  =  DBHT 

VAR(12)  =  DENT 

VAR(13)  =  HTST 

VAR(14)  =  BAST 

JDENT     =  DENT 

5.  Replace  the  statement  labeled  55  with  two 
statements: 

55  IF(K  .EQ.  1)  GO  TO  56 
CALL  PONCUT 

This  will  bypass  the  initial  thinning. 

6.  Place  three  statements  between  GO  TO  60 
and  the  statement  labeled  58: 

56  ADDHT  =  HTST  -  HTSO 
HTCUM  =  HTCUM  +  ADDHT 
GO  TO  59 

7.  Add  a  label  to  the  statement  for  STAND, 
three  statements  after  label  58: 

59  STAND  =  DENT 


8.  Place  four  statements  just  before  the  CON- 
TINUE statement  labeled  160: 
DBHT  =  VAR(ll) 
DENT  =  VAR(12) 
HTST  =  VAR(13) 
BAST  =  VAR(14) 
JDENT  =  VAR(12) 


To  Add  Variability: 

A  modification  of  PONYLD  primarily  of 
interest  as  a  research  tool  is  the  addition  of 
random  elements  to  each  prediction.  Several 
stands  with  exactly  the  same  basal  area,  site 
index,  etc.  do  not  all  produce  the  same  periodic 
growth  or  have  the  same  volume  conversion 
factors.  A  prediction  equation  will,  however, 
indicate  the  same  average  response  for  all  the 
stands.  Response  in  the  presence  of  variability 
can  be  studied  by  adding  statements  to  compute 
the  random  elements,  running  the  program 
several  times,  and  then  analyzing  the  results 
statistically. 

Modification  of  PONYLD  to  add  variability 
to  computations  of  DBHO  illustrates  how  such 
changes  can  be  made.  The  following  statements 
are  added  to  the  program  immediately  after  the 
statements  that  compute  DBHO  and  round  it 
to  0.1  inch. 

98    IDIV  =  (17.0  *  GNTR  +  3.0)  /  1024.0 
NGNTR  =  GNTR 

GNTR  =    (17  *   NGNTR  +  3)  -  1024 

*  IDIV 

IF  (GNTR  .GT.  1000.0)  GO  TO  98 
IF(GNTR  .LT.  0.0)  GO  TO  98 
Al  =  GNTR  /  100.00 
A2  =  Al  *  Al 

RES  =  0.9565  *  Al  -  0.0523  *  A2  -  0.0063 

*  Al  *  A2  +  0.00084  *  A2  *  A2 
-  3.3009 

IRES  =  RES 

IF(RES  .LT.  0.0)  IRES  =  RES  -  0.5 
IF(RES  .GT.  0.0)  IRES  =  RES  +  0.5 
ADJ  =  IRES 

DBHO  =  DBHO  +  ADJ  *  0.1 


The  statement  for  IDIV  is  a  congruential 
pseudo-random  number  generator  (Greenberger 
1961).  A  value  of  the  variable  GNTR  is  read 
into  memory  at  the  beginning  of  the  program. 
It  is  the  "seed"  of  the  generator  and  may  be 
any  number  from  0  to  1023,  inclusive.  This 
variable,  or  the  entire  statement,  is  varied  in 
successive  runs  of  PONYLD  to  give  the  results 
an  opportunity  to  differ. 
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The  statement  for  RES  is  an  empirical  distri- 
bution function  obtained  by  fitting  a  polynomial 
to  the  normally  distributed  residuals  of  the 
DBHO  equation  (Evans  et  al.  1967). 

The  computation  is  made  in  three  steps: 

(1)  generation  of  a  pseudo-random  number, 

(2)  use  of  this  number  as  an  independent 
variable  to  compute  the  value  of  a  normally 
distributed  residual  (range:  -0.3  to  +0.3  inch 
for  DBHO),  and  (3)  addition  of  the  residual 
algebraically  to  the  computed  value  of  DBHO. 

To  Include  Overstory: 

An  addition  may  be  needed  for  tables  with 
shelterwood  or  seed  tree  cuts  if  the  regeneration 
period  is  long  enough  for  the  next  crop  to 
affect  growth  of  the  overstory.  In  this  case, 
PONYLD  should  be  run  twice.  The  first  run 
provides  basal  areas  at  appropriate  times  early 
in  the  life  of  a  stand.  Changes  are  then  made 
for  the  second  run,  which  will  be  the  pro- 
duction run.  BAST  of  the  DBHO  equation 
becomes  total  basal  area  of  overstory  and  under- 
story  for  growth  projections  made  late  in  the 
life  of  the  overstory  after  the  new  crop  trees 
are  large  enough  to  influence  overstory  growth. 

Description  of  Program  PONCHK 

Most  functions  described  in  the  section 
headed  Information  Needed  are  computed  with 


rather  large  amounts  of  plot  data.  The  func- 
tions must  apply  to  many  combinations  of 
several  stand  variables  and  to  the  wide  range 
of  stand  density  possible  in  thinned  and  un- 
thinned  stands.  It  is  possible,  however,  to 
calculate  additional  data  for  use  in  regression 
analyses  of  several  of  these  functions  from 
lesser  amounts  of  actual  field  data.  The  de- 
rived functions  must,  of  course,  be  checked 
for  accuracy.  If  acceptable,  they  may  be  used 
in  production  runs  of  PONYLD.  If  not  good 
enough  for  this,  they  may  be  used,  until  addi- 
tional field  data  are  obtained,  for  preliminary 
modeling  and  adapting  PONYLD  to  local  equip- 
ment and  species. 

Program  PONCHK,  listed  in  appendix  3, 
uses  plot  stand  tables  to  calculate  volumes, 
volume  conversion  factors,  and  changes  in 
diameter  and  height  due  to  thinning.  Thin- 
nings that  leave  about  10,  17,  20,  25,  33,  50, 
67,  and  75  percent  of  the  trees  are  simulated. 
Printed  output  of  PONCHK  provides  both 
dependent  and  independent  variables  needed 
to  obtain  functions  for  CUFT,  FCTR,  PROD, 
DBHE,  and  ADDHT.  Also  shown  are  the 
volumes  removed  by  each  level  of  thinning, 
including  the  merchantable  cubic  volume  in 
trees  below  saw  log  size. 

Data  cards  needed  to  run  PONCHK  are 
shown  in  the  accompanying  tabulation.  The 
last  card  of  the  data  deck  is  a  type  1  card 


Order  and  Contents  of  the  Data  Deck  for  Program  PONCHK 

Card      Number  of  Variable 
type  cards         name         Columns       Format       Description  of  variable 


1  per 
plot 


1  per 
tree 


NPLT 
NTIM 

NWHN 

AREA 
AGE 

DBH 
HT 

CC 


1-3 
4-5 

6-  7 

8-12 
13-16 

1-3 
4-6 

7-  9 


13 
12 

12 

F5.3 
F4.0 

F3.1 
F3.0 

F3.0 


Plot  number 

Diameters  from  measured 
d.b.h.  (2)  or  from  incre- 
ment borings  (1) 
Field  plot  is  unthinned 
(1)  or  thinned  (2) 
Plot  area 

Average  age  of  main 
stand 

Tree  diameter  to  0.1  inch 
Tree  height  to  the  near- 
est foot 

Tree  crown  class:  domi- 
nant (1),  codominant  (2), 
intermediate  (3),  or  over- 
topped (4) 


c 
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with  888  punched  in  the  NPLT  field.  Type  2 
cards  of  each  plot  are  arranged  in  order  of 
increasing  diameter.  The  last  type  2  card 
of  each  plot  has  99.9  (punched  999)  in  the 
DBH  field. 

The  sequence  of  operations  is  shown  by 
the  COMMENT  statements  of  the  program. 
Summations  and  computations  of  averages 
follow  usual  procedures.  Simulation  of  thin- 
nings by  PONCHK  is  accomplished  for  a  plot 
as  follows: 

1.  A  modulo  for  the  pseudo-random  number 
generator,  which  is  appropriate  for  the  num- 
ber of  trees  per  acre,  is  computed. 

2.  As  many  pseudo-random  numbers  are  gener- 
ated as  there  are  trees  on  one  acre,  using 
the  generator  provided  and  the  "seed"  value 
of  the  statement  labeled  1  (Greenberger  1961). 

3.  Trees  are  rearranged  at  random. 

4.  Groups  of  randomly  arranged  trees  are  created 
with  the  size  of  each  group  determined  by 
the  percentage  of  trees  to  be  retained. 

5.  Appropriate  trees  from  each  group  are  tallied. 
For  example,  with  PRET  =  25.0,  the  largest 
tree  of  each  group  of  four  is  tallied.  For 
PRET  =  75.0,  the  smallest  tree  of  each  group 
of  four  is  ignored  and  the  other  three  are 
tallied. 

6.  Volumes  and  other  values  are  computed 
for  each  percentage  of  trees  retained. 

In  tests  to  date,  functions  computed  from 
data  from  simulated  thinnings  agree  with  results 
of  field  operations.  For  example,  equations  for 
DE-HE  derived  from  simulated  data  usually 
predict  post-thinning  diameters  that  vary  0.1 
inch  or  less  from  actual  values. 
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APPENDIX  1 
Listing  of  Program  PONYLD 


C   TO  COMPUTE   AND  PRINT   YIELD  TABLES   FOR   MANAGEO  EVEN- AGED  STANDS. 
C 

C   DEFINITIONS    OF  VARIABLES. 

C 

C  ADDHT  =   INCREASE   OR  DECREASE   IN  AVERAGE  STAND  HEIGHT   BY  THINNING. 

C  AGEO  =    INITIAL    AGE    IN   YIELO  TABLE. 

C  BASC  =  BASAL   AREA  CUT   PER  ACRE. 

C  BASO  =   BASAL    AREA   PER    ACRE    BEFORE  THINNING. 

C  BAST  =  BASAL   AREA  PER   ACRE   AFTER  THINNING. 

BOARD  FEET  CUT   PER  ACRE. 
BOARD   FEET   PER   ACRE   BEFORE  THINNING. 
BOARD  FEET   PER   ACRE   AFTER  THINNING. 
MERCHANTABLE  CU.   FT.   CUT   PER  ACRE. 
MERCH.   CU.    FT.    PER   ACRE  BEFORE  THINNING. 
MERCH.   CU.   FT.    PER  ACRE   AFTER  THINNING. 
COMBF    =    MINIMUM  COMMERCIAL   CUT,    BOARD  FEET. 
COMCU  =  MINIMUM  COMMERCIAL  CUT,    CU.  FT. 
DBHO  =    AVERAGE    STAND  O.B.H.    BEFORE  THINNING. 

AVERAGE    STAND   D.B.H.    AFTER  THINNING. 
TREES  CUT   PER  ACRE. 
TREES   PER   ACRE   BEFORE  THINNING. 
TREES   PER  ACRE   AFTER  THINNING. 

PERCENTAGE,    AS    A  DECIMAL,    OF   TREES    THAT   DIE   DURING  PERIOD 
R I  NT. 

GROWING  STOCK  LEVEL  FOR   INTERMEDIATE  CUTS  AFTER  FIRST. 
DSTY  =  LOWEST   VALUE  OF  DLEV  USEO   IN  A  TEST. 

G IDE   =  BASE  FOR  GROWING  STOCK   LEVELS,    80.0   IN  EXAMPLE  SHOWN. 
HTSO  =  TREE  HEIGHT  BEFORE  THINNING. 
HTST  =  TREE  HEIGHT  AFTER  THINNING. 

INVL(l)    =   NEW   CUTTING   CYCLE   AFTER    REGENERATION  CUT  [. 
JCYCL  =   INTERVAL  BETWEEN   INTERMEDIATE  CUTS. 

JSBD   =    SUM  OF    BOARD  FEET    FROM   ALL   CUTS   WITH  YIELD  OF   COMBF  OR 
LARGER. 

JSMC  =   SUM  OF   MERCH.   CU.   FT.   FROM   ALL  CUTS  WITH  YIELD  OF  COMCU  OR 
LARGER. 

JSTF  =   SUM  OF   TOTAL  CU.   FT.    FROM  ALL  CUTS. 
MIX  =  NUMBER  OF   STOCKING  LEVELS  EXAMINED   PER  TEST . 
NTSTS  =  NUMBER  OF   TESTS   PER  BATCH. 

PRET  =  PERCENTAGE   OF  TREES  RETAINED  AFTER  THINNING. 
REGNII)    =   STANO  AGE  WHEN  REGENERATION  CUT   I  OCCURS. 
RI NT   =  NUMBER   OF    YEARS    FOR  WHICH   PROJECTION    IS  MADE. 
FINAL    AGE    IN  YIELD  TABLE. 
SITE  INDEX. 

GROWING   STOCK  LEVEL  FOR    INITIAL  THINNING. 
TOTAL  CUBIC   FEET  CUT  PER  ACRE. 
TOTAL  CUBIC   FEET  PER   ACRE   BEFORE  THINNING. 
TOTAL  CUBIC  FEET   PER  ACRE   AFTER  THINNING. 

PERCENT  OF  PREVIOUS  DLEV  TO  BE  LEFT  AT  REGNII),  ENTERED 
DECIMAL. 


BDFC 
BDFO 
80FT 
CFMC 
CFMO 
CFMT 


DBHT 
OENC 
DENO 
OENT 
DIED 

DLEV 


ROTA 
SITE 
THIN 
TOTC 
TOTO 
TOTT 
VLLVt I ) 


COMMON  BA, BAST, CUFT, DBHO, DBHT, DENO, FCTR,HITE.GIDE,PRET,PROD,REST,S 
1 T AND , VDM 

DIMENSION   INVL ( 3)  * REGN ( 3 ) ,VAR( 10) , VLLVt 3) 


AL0G101 SI TE ) 


C   PROVIDE   FOR   SEVERAL   GROWING  STOCK   LEVELS  PER  TEST. 

C 

30  DO   160  M=1,MIX 
ADDHT   =  0.0 
BOFO   =  0.0 
BDFT  =  0.0 
CFMO  =  0.0 
CFMT  =  0.0 
HTCUM   =  0.0 
JSBD  =  0 
JSMC   =  0 
JSTF  =  0 
TEM   =  M 
DLEV  =  (DSTY 
BASO  =  DENO  * 
C 

C   OBTAIN   AVERAGE   HEIGHT    ANO   VOLUMES    PER  ACRE. 

C 

C  STATEMENTS  FOR  HTSO  AND  IF  STATEMENT   ARE   SPECIES- SPECIF IC. 

C 

IF ( AGEO   -GT.   55.01   GO  TO  35 

HTSO  =  0.01441   *  AGEO  *  SITE  -  0.12162  *  AGEO  -  1.50153 
GO  TO  40 

35  HTSO  =  0.59947  -  61.5019  /   AGEO  +  0.80522 
18  *  ALOGIO(SITE)    /  AGEO 

HTSO    =    10.0   **  HTSO 
<.0  HITE  =  HTSO 

OA   =  BASO 

STAND   =  DENO 

VDM  =  DBHO 

CALL  PONVOL 

TOTO  =  CUFT 

BDFO  =  CUFT  *  PROD 

CFMO  =  CUFT  *  FCTR 

REST  =  THIN 

C 

C   ENTER  LOOP  FOR  REMAINING  COMPUTATIONS   AND  PRINTOUT. 

C 

DO   130  K  =  l  ,100 

C 

C   CHANGE   STANDAROS    IF   A  REGENERATION  CUT   IS  OUE. 
C 

43   IF ( AGEO   .GE.   ROTA)   GO  TO  60 

1FIAGE0   .LT.   REGN ( 1 ) )   GO  TO  55 
I F ( AGEO   . NE.    REGN (1)1   GO  TO  45 
DLEV   =   DLEV   *  VLLV(l) 
REST  =  DLEV 
JCYCL    =    INVL ( 1 ) 
GO  TO  55 

45   I F ( AGEO    .NE.    REGN ( 2  1  )   GO  TO  50 
DLEV  =  DLEV  *  VLLV12) 
REST  =  OLEV 
JCYCL   =  INVL(2) 
GO  TO  55 

50   I F I AGEO   .NE.    REGN(3)>   GO  TO  55 
DLEV  =  DLEV  *  VLLVI3) 
REST  =  OLEV 
JCYCL    =  INVL13I 

C 

C    INCREASE   D.B.H.    BY   THINNING   AND  COMPUTE    POST-THINNING  VALUES. 


REAO    15,5)   NTSTS, GIDE 
5   FORMAT    ( I4.F4.0) 

1F(NTSTS  .LE.  0)  GO  TO  170 
IF1GIDE    .LE.   0.0)    GO  TO  170 


READ    (5,10)    COMCU, COMBF 
10   FORMAT  (10F8.3) 
VARtB)    =  COHCU 
VAR I  9 1    =  COMBF 


C    EXECUTE    PROGRAM  ONCE    FOR   EACH   SET    OF    INITIAL    VALUES   OF  INTEREST. 


55  CALL  PONCUT 

JDENT  =  (BAST  /  (0.0054542  *  DBHT  «  DBHT 11 
DENT  =  JDENT 

BAST  =  0.0054542  *  DBHT  *  DBHT  *  OENT 

I F I  BAST   .LT.    BASO)   GO  TO  58 

BAST  =  BASO 

HTST  =  HTSO 

OENT  =  DENO 

JDENT   ~    DENO  ♦  0.5 

DBHT   =  DBHO 

TOTT  =  TOTO 

BDFT  =  BDFO 

CFMT  =  CFMO 

GO  TO  60 

C 

C  STATEMENT   FOR   ADDHT   IS   SPEC  IE S- SPECIFIC. 


C    REAO  CUTTING    INTERVAL    AND   LEVELS  PER    TEST   FROM  CARD   TYPE  THREE. 

C 

READ    (5,15)    JCYCL, MIX 
15   FORMAT  (214) 

IF1JCYCL   .LE.   0   .OR.   MIX   .LE.   0)    GO  TO  170 

J  TEM  =  JCYCL 

C 

C   READ    INITIAL    STAND  VALUES    FROM  CARD   TYPE  FOUR. 

C 

READ    15,10)    AGEO, DBHO, DENO, DSTY, RINT, SITE, THIN 

VAR(l)    =  AGEO 

VAR ( 2 1   =  OBHO 

VAR ( 3 )   =  DENO 

VARI4)    =  DSTY 

VAR<5)    =  RINT 

VARI6)    =  SITE 

VARI7)    =  THIN 

C 

C    READ   SILVICULTURAL    CONTROLS   FROM  CARD   TYPE  FIVE. 

C 

READ    (5, 10) REGN  1 1  )  ,VLLV( I  ) ,  INVLl  1 ) , REGN ( 2  ) , VLL V ( 2  ) , INVL 12  ) , REGN I  3 ) 
VAR(  10)    =  REGN ( 1 ) 
00  20  L=l , 10 

IF(VARIL)    .LE.   0.01   GO  TO  170 
20  CONTINUE 

DLEV  =  0.0 
DO  25  NA=l,3 
L   =  4  -  NA 

IFIREGN(L)    .EQ.   0.0)   GO  TO  25 
ROTA  =  REGN(L) 
GO  TO  30 
25  CONTINUE 


58  AODHT  =   7.64833  -  3.82286   *  ALOGIO(PRET) 
HTCUM  =  HTCUM  +  ADDHT 
HTST  =   HTSO  +  ADDHT 
STAND  =  DENT 
VDM  =  DBHT 
HITE   =  HTST 
BA  =  BAST 
CALL  PONVOL 
TOTT  =  CUFT 
BDFT  =  CUFT  *  PROD 
CFMT   =  CUFT  *  FCTR 

C 

C   CHANGE    MODE    AND  ROUND   OFF    FOR  PRINTING. 


60  JAGEO 
JSITE 
JDENO 
JHTSO 
JTOTO 
JTOTO 
JBASO 
JCFMO 
JCFMO 
JBDFO 
JBOFO 
JHTST 
J  TOTT 
JTOTT 
JCFMT 
JCFMT 


AGEO 
SITE 

DENO   ♦  0.5 
HTSO  +  0.5 
(TOTO    *  0.11    +  0.5 
JTOTO  *  10 
BASO  ♦  0.5 
(CFMO  *  0.1)    +  0.5 
JCFMO    *  10 
I  BDFO    *   0.01)    ♦  0.5 
JBDFO   *  100 
HTST   ♦  0.5 
(TOTT  *  0.1) 


0.5 


JTOTT 
I  CFMT 
JCFMT 
IFIJCFMT  .GT. 
JBDFT  =  (BDFT 
JBDFT  =  JBDFT 
IFUHDFT  .GT. 
J0AST  =  BAST 
JDENC   =   JDENO  -  JDENT 


10 

O.I)    +  0.5 
10 

JCFMO)  JCFMO  = 
»  0.011  +  0.5 
*  100 

JBDFO)  JBOFO  = 
■  0.5 
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JBASC  =  JBASO  -  JBAST 
JTOTC  =  JTOTO  -  JTOTT 
JCFMC  =  JCFMO  -  JCFMT 
IFIJCFHC  .LE.  0)  JCFMC  *  0 
JflDFC  =  JBDFO  -  JBOFT 
IFIJBDFC   .  LE  .   0)    JBOFC   ■  0 

C 

C   SUM  PERIODIC   CUTS  FOR   LAST  LINE  OF  YIELD  TABLE. 

C 

IFUGEO   .  GE .    ROTA)   GO  TO  70 
JSTF   =  JSTF  ♦  JTOTC 
CFMC   =  JCFMC 

IFICFHC   . LT.   COMCU)   GO  TO  65 
JSMC   =   JSMC  ♦  JCFMC 
65  BDFC   =  JBDFC 

I F ( BDFC    .LT.    COHBF)    GO  TO  70 
JSBO  =   JSBD  ♦  JBDFC 

C 

C  WRITE  HEADINGS  FOR   YIELD  TABLE. 
C 

70  IFIK   ,GE.   2)   GO  TO  92 

C 

C  CHANGE   TABLE   HEADING  FOR  OTHER  SPECIES. 

C 

WRITE    16,80)  JSITE.JCYCL 
80  FORMAT    (  1 H 1  , / // , 3 7X , 6 2HY I  ELDS  PER   ACRE  OF  MANAGED,   EVEN- AGEO  STANO 
IS  OF   PONDEROS A  PI NE/ IH0.48X ,  1 1HS I TE   INDEX   , 13, IH, , 14.  19H-YEAR  CUTT 
2ING  CYCLE  I 
WRITE    (6,82)    TH I N, DLE V 
82  FORMAT    t 1H0.41X.26HTHINNING  LEVEL S=    INITIAL  -.F6.0.14H,  SUBSEQUENT 
1  -,F6.0) 
WRITE  (6,84) 

84  FORMAT    ( IH0.25X.38HENTIRE   STAND  BF  FORE   AND  AFTER  TH INN  I  NO, 28X , 26HP 
1ERI0DIC    INTERMEDIATE  CUTS) 
WRITE  (6,86) 

86   FORMAT    ( I HO , 9X , 5H STAND , 10X , 5HBAS AL , 3X , 7H AVE RAGE, 2X , 7HAVER AGE , 3X , 5H 
1TOTAL,3X,9HM£RCHANT-,3X,9HSAWTIM(JER,9X,5HBASAL,4X,5HTOTAL,3X,9HMER 
2CHANT-,3X,9HSAWTIHBER) 
WRITE  (6,88) 

88  FORMAT    1 1H    , 1 0 X , 3H AGE ,4X , 5HT REES , 3 X , 4 H ARE A, *X , 6HD . B  . H . , 3X , 6HHE I GH T 
1,2X,6H  VOLUME,  2  X.UHABLE  VOLUME,  4  X  ,  6HV0LUME,  3X  ,  5HTREE  S»  3X,  4HAREA  ,  3X 
2.6HV0LUME.2X, 1 1HABLE   VOLUME »4X ,6H VOLUME ) 
WRITE  (6,90) 

90  FORMAT  (1H  , 8 X , 7 H ( YE ARS ) , 3X , 3HN0 . , 3X , 6HS0 . F T . , 4X , 3H IN . , 6X , 3HF T . , 4 X 
lt6HCU.FT.,5X,6HCU.FT.t6X.6HBD.FT.,4X,3HN0.,3X,6HSG.FT.,2X,6HCU.FT. 
2,5X,6HCU.FT. ,6X,6HBD.FT. ) 

C 

C   WRITE   TABLE   ENTRIES  OF   DIAMETER,   VOLUMES,  ETC. 
C 

92   WRITE    (6,94)    JAGEO, JDENO, JB ASO, 08H0, J HTSO , J TOTO , JCFMO , JBDFO 
94  FORMAT    (  IHO , 9X , I  4 , 4X , I5,2X,14,5X,F5.1,5X,I3,4X,15,6X,I5,6X,I6) 
I F 1 AGEO   .GE.    ROTA)   GO  TO  135 

WRITE    16,96)    JAGEO, JOENT , JBAST , DBHT , JHTST , JTOTT, JCFMT , JBDF  T , JDENC , 
LJflASC, JTOTC, JCFMC, JBOFC 
96  FORMAT    (LH    , 9 X , 1 4 , 4 X , I  5 , 2 X , I  4 , 5 X , F5 .  1 , 5X ,  I  3 , 4X ,  I  5 , 6X , I  5 , 6X  ,  I  6  ,  4X ,  I 
l5,3X,I3,5X,I4,6X,I4,8XtI5) 

C 

C  COMPUTE   VALUES  FOR  EACH  PERIOD.    THIN   AS  SPECIFIED. 

C 

IR  INT   =   R  INT 

IK  =  JCYCL   /  IRINT 

DO   120  L=1,IK 

AGEO  =  AGEO  +  RINT 

I F ( AGEO  .GT.   ROTA  I   GO  TO  135 

C 

C  COMPUTE  NEW  D.B.H.   BEFORE  THINNING  AND  ROUND  OFF  TO  0.1  INCH. 
C 

C  STATEMENT   FOR  DBHO   IS   SPEC IE S~ SPECIFIC. 

C 

DBHO  =   1.0097  *  DBHT  +  0.0096  *  S I T E- ( 1 . 5766* ALOG 10 < B A S T ) ) + 3 . 3021 
I DBHO  =   OBHO  *   10.0  ♦  0-5 
OBHO  =   I DBHO 
DBHO   =   DBHO   *  0.1 

C 

C  STATEMENT  FOR  DIED   IS   SPECIES- SPECIFIC.   CHANGE   10.0   IN  IF 

C  STATEMENT   IF  DIED  STATEMENT   APPLIES  TO  LARGER  TREES. 

C 

IF108HT    .GE .    10.0)   GO  TO  100 

DIED  =  0.00247  +  0.00124  *  DBHT   +   0.00028  *  OBHT  *  OBHT   ♦  0.000005 
121  »  BAST  *  BAST  -  0.0000905  *  DBHT  *  BAST 
IFIDIED   .LT.    0.0)   OIEO  =  0.0 
DENO  =  DENT   *    (1.0  -  DI60I 
MNK  =  OENO  ♦  0.5 
OENO  =  MNK 
GO  TO  105 
100  DENO  =  DENT 

105  BASO  =  DENO  *    (0.0054542   *  DBHO  *  DBHO) 

C 

C  OBTAIN  AVERAGE  HEIGHT   AND  VOLUMES  PER  ACRE. 
C 

C  STATEMENTS  FOR  HTSO  AND  IF  STATEMENT   ARE   SPECIES- SPECIFIC. 

C 

I F I AGEO   .GT.    55.0)   GO  TO  110 

HTSO  =   0.01441   *  AGEO  *  SITE  -  0.12162  *   AGEO  -  1.50953 
GO  TO  115 

110  HTSO  =  0.59947  -  61.5019  /  AGEO  +  0.80522  »  ALOGIO(SITE)  ♦  20.5252 
18  *  ALOGIO(SITE)    /  AGEO 

HTSO  =  10.0  **  HTSO 
115  HTSO  =  HTSO  ♦  HTCUH 

STAND  =  DENO 

VOM  =  DBHO 

HITE   =  HTSO 

BA  =  BASO 

CALL  PONVOL 

TOTO  =  CUFT 

BDFO  =  CUFT   *  PROD 

CFMO  =  CUFT   6  FCTR 

C 

C   TEST   IF  REGENERATION  CUT   IS  DUE. 
C 

DO    118  KU=1,3 

I F ( AGEO   .EO.   REGNIKU))   GO  TO  43 
118  CONTINUE 

C 

C   CHANGE   MODE    AND  ROUND  OFF    FOR  PRINTING. 
C 

IFtl   . EO.    IK)   GO  TO  125 

KDENO  =  DENO  +  0.5 

K AGEO  =  AGEO 

KHTSO   =   HTSO    ♦  0.5 

KBASO  =   BASO  +  0.5 

KTOTO  =    (TOTO  *  0.1)   ♦  0.5 

KTOTO  =   KTOTO  *  10 

KCFMO  =    (CFMO  *  0.1)    ♦  0.5 

KCFMO  =   KCFMO  *  10 

KBOFO  =    ( BDFO  *  O.Ol)    +  0.5 

KBDFO  =   KBDFO   *  100 


C 

C   WRITE  VALUES  FOR   THE   PERIOD   IF  THINNING   IS  NOT  DUE. 

C 

WRI TE    (6,94)    K AGEO, KDENO, KB ASO, DBHO, KHTSO, KTOTO, KCFMO, KBDFO 

DBHT    =  DBHO 

BAST  =  BASO 

DENT  =  DENO 
120  CONTINUE 
125  REST  =  DLE V 
130  CONTINUE 

C 

C  ADD  FINAL  CUTS    TO    TOTAL   YIELDS    AND  WRITE    TOTAL  YIELOS. 

C 

135  JSTF  =   JSTF   ♦  JTOTO 
CFMO  =  JCFMO 

IF1CFM0  .LT.   COMCU)   GO  TO  140 

JSMC  =  JSMC  «•  JCFMO 
140  BOFO  =  J80F0 

I F (BDFO  .LT.   COMBF)   GO  TO  145 

JSBO  =  JSBD  +  JBDFO 
145  WRITE    (6,150)    JSTF , JSMC , JSBD 

150  FORMAT    (IHO, /,67X,12H TOTAL   Y I  EL OS , 20X ,  1 4 , 6X , 1 4 , 8X ,  I  5  ) 
WRITE    (6,155)   COMCU, COMBF 

155  FORMAT    (  IHO, // , 1  IX ,44HMI N I  HUM  CUTS  FOR    INCLUSION    IN  TOTAL    YIELDS — 
I.F6.0.15H  CUBIC  FEET   AND,F7.0,11H  80AR0  FEET) 

WRITE  (6,156) 

156  FORMAT   ( IHO , 1  OX , 66HMERCH.   CU.   FT.   -  TREES  6.0   INCHES  D.B.H.   ANO  LA 
1RGER    TO  4-INCH   TOP. ) 

WRITE  16,157) 

157  FORMAT    (  IHO,  10X , 60HBD .    FT.    -   TREES    10.0    INCHES   D.B.H.    AND  LARGER  T 
10   8-INCH    TOP. ) 

C 

C   PREPARE   FOR  NEXT   TABLE   OF   THE  TEST. 

C 

AGEO  =  VAR( 1 ) 

DBHO   =  VAR(2) 

DENO  =  VARI3) 

JCYCL  =  JTEM 
160  CONTINUE 

GO  TU  200 
170  WRITE  (6,175) 

175  FORMAT    ( 1H 1 , / / / , 1 OX , 66HE XECUT I  ON  STOPPED   BECAUSE  OF  NEGATIVE  OR  ZE 

1R0    ITEM  ON  A  DATA  CARD.) 
200  CALL  EXIT 
END 


SUBROUTINE  PONCUT 

C 

C   TO  ESTIMATE    INCREASE    IN  AVERAGE  D.B.H.   DUE  TO  THINNING. 

C 

COMMON  BA, BAST, CUFT, OBHO, DBHT, DENO, FCTR, HITE.GI 06, PRET, PROD, REST, S 
1TAND.VDH 

C 

IF(08H0    .LT.    9.4)    GO   TO  30 

C 

C   COMPUTE    D.B.H.    IF   DBHO    IS   LARGE    ENOUGH   FOR   BASAL    AREA    TO  REMAIN  CONSTANT. 

C 

PRE  T  =  1C0.0 
00  21  KJ=i,100 

C 

C  STATEMENTS  FOR   DBHE    AND  PDBHE    ARE  SPECIES-SPECIFIC. 

C 

I F ( PRE T   .LT.   50.0)   GO   TO  5 

OBHE   =  0-73365  +   1.02008   *  DBHO  -  0.01107  *   ( PRET  -  50.0)   -  0.0001 
14  *    i PRE T  -  50.0)    *   (PRET  -  50.0) 
GO  TO  11 

5  PDBHE   =   0.49401   +  0.71890  *   AL0G1 0 1 DBHO  1   -  0.22530  *  AL0G101 PRET ) 
1   +  0.12616  *  AL0G1 0 ( DBHO)    *  AL0G10 ( PRET ) 
DBHE   =   10.0  *»  PDBHE 
11   IOBHE   =  OBHE     *    10.0  +  0.5 
DBHE   =   I DBHE 
DBHE    =  DBHE   *  0.1 
DENE   =  DENO  *  PRET  *  0.01 
NOENE  =  OENE   ♦  0.5 
DENE   =  NOENE 

BASE    =    0.0054542    *   DBHE   *   DBHE   *  DENE 

NBASE   =   BASE   *  10.0  +  0-5 

BASE   =  NBASE 

BASE   =   BASE  *  0. 1 

TMPY  =  0.0054542   *   DBHE  *  DBHE 

TEM  =  BASE  -  REST 

IFITEM   .LE.    TMPY)    GO  TO  70 

IFITEM    .LT.    4.0)    GO  TO  20 

PRET  =  PRET  -  1.0 

GO  TO  21 

20  PRET  =   PRET  -  0.3 

21  CONTINUE 
GO   TO  70 

C 

C  COMPUTE  O.B.H.    IF   BASAL  AREA   INCREASES  WITH  D.B.H. 
C 

30  PRET  =  40.0 

I F ( DBHO    .GT.    7.0)    PRET   =  70.0 
DO  65  J=l,100 

C 

C  STATEMENTS  FOR  OBHE  AND  PDBHE  ARE  SPECIES-SPECIFIC. 

C 

I F ( PRET  . GE.   50.0)   GO  TO  40 

PDBHE    =   0.49401    +   0.71890  *    AL0G1 0 (DBHO )    -   0.22530   *  ALOGIO(PRET) 


1   +   0.12616  *   AL0G1 0 1 DBHO)    *   AL0G1 0  I PRET ) 
DBHE    =    10.0   **  PDBHE 
GO  TO  45 

40  OBHE   =  0.73365  *   1.02008   *  DBHO  -  0.01107  *   (PRET  -  50.0)   -  0.0001 

14  *    I PRET  -  50.0)    *   I PRET  -  50.0) 
45  I DBHE  =  OBHE   *   10.0  +  0.5 

DBHE  =   I DBHE 

DBHE  =  DBHE   *  0.1 

DENE   =  OENO  *    (PRET  *  0.01) 

NDENE   =  DENE   +  0.5 

OENE  =  NOENE 

BASE  =  0.0054542   *  DBHE  *  DBHE   *  DENE 
NBASE   =   BASE   *   10.0  ♦  0-5 
BASE   =  NBASE 
BASE   =   BASE   *  0.1 

C 

C  CHANGE   STATEMENTS  FOR  BREAK,    BUST   ANO  FIRST   3  STATEMENTS  FOR  DBHP 

C  IF   OTHER  GROWING   STOCK  LEVEL   BASE   THAN   TABLE    1   IS  USED. 

C 
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BREAK   =  49.9    *  REST   /  GIOE 
I F ( BASE    .GT.    BREAK)   GO  TO  50 
OBHP   =    IGIDE   /  REST )    *    (0.08682  * 
GO  TO  52 

50  BUST  =  66.2  *    I  RE  ST   /  GIOE) 
I F ( BASE    • GT.    BUST)   GO  TO  51 

OBHP  =  IGIDE  /  REST)  *  (0.10938  * 
GO  TO  52 

51  TMPY  =   BASE  *    (GIOE   /  REST) 
TEH  =  THPY  *  THPV 

OBHP   =    19.04740   *  TMPY  -  0.26673 
I  -  448.76833 
I F ( THPY   .GT.   GIOE)    OBHP  =   OBHO  ♦ 

52  IDBHP  =  OBHP   *   10.0  ♦  0.5 
OBHP  =  IDBHP 

DBHP  =   OBHP  *  0.1 
IFIDBHP  -   DBHE)  60,70,61 

60  PRET  =   PRET  *  1.02 
GO  TO  65 

61  PRET   =   PRET  *  0.9B 
65  CONTINUE 

70  DBHT  =  D8HE 

C 

C   COMPUTE    POST-THINNING    BASAL  AREA. 


C   COMPUTE    TOTAL  CUBIC   FEET    PER  ACRE. 

C 

D2H  =  VDM  *   VOM  *  HITE 


 STATEMENTS  FOR  CUFT  ANO   IF   STATEMENT   ARE  SPECIES-SPECIFIC. 

IFI02H   .GT.   6000.0)   GO  TO  5 

CUFT   =    (0.00225   *   D2H  -   0.00074    *    BA   *  0.03711)    *  STAND 
GO  TO  10 

5  CUFT  =    (0.00247  *  D2H  +  0.00130  *   BA  -   1.40286)   *  STAND 
10    IFIVDM    . LT .    5.0)    GO  TO  40 


VOLUMES  TO  4.0-INCH  TOP 


I F (DBHT   .GT.    5.0)   GO  TO  75 

SOFT  =    U. 58495   *  DBHT  -  11.09724 

GO  TO  76 

75  IF10BHT   . GE.    10-0)   GO  TO  77 
TEH  =  DBHT  *  DBHT 

SOFT  =  7.76226  *  DBHT   +0.85289  *  TEM 

76  BAST  =    (REST   /  GIOE)    *  SOFT 
GO  TO  80 

77  BAST  =  REST 
80  RETURN 

END 


DBHT-3. 45624 


IFIVDM   .GT.    6.7)    GO  TO  15 
FCTR  =  0.26612  *  VOH  -  1.12689 
GO  TO  25 
15   IFIVDM  .GT.   10.4)   GO  TO  20 

FCTR  =  3.46993  -  0.12017  *  VDH  -  13.41984  /  VDM 
GO  TO  25 

20  FCTR  =   0.99666  -   0.66932    /  VDH 
25    IFIVOH    .LT.    8.0)    GO  TO  40 


SUBROUTINE  PONVOL 


C    TO  COHPUTE    VOLUHES   PER    ACRE    IN  VARIOUS  UNITS. 


C  OBTAIN  CONVERSION  FACTORS  FOR  BD.  FT.  -  VOLUMES  TO  B-INCH  TOP  IN  TREES 
C   10.0    INCHES  D.B.H.    ANO  LARGER. 

C 

C  STATEMENTS  FOR   PROO  AND  IF  STATEMENT   ARE  SPECIES-SPECIFIC. 


IFIVDM   .GT.    11.9)'  GO  TO  30 
PROD  =  0-87783   *  VDH  +  0.00660 
GO  TO  40 

30  PROD  =   5.10752  +  0.10712   *  VDH 
40  RETURN 
END 


BA  -  7.27957 

0.00185  *  BA  -   36.20229  /  VDM 
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APPENMX  1 
Output  of  PONYLD 


S  fANO 
AGE 
(YEARS) 

30 
30 


50 
50 


70 
70 


950 
288 


284 
172 


171 
104 


YIELDS   PER   ACRE  OF   MANAGEO,    EVEN-AGEO   STANDS  OF   PONDEROSA  PINE 
SITE   INDEX      70,      20-YEAH  CUTTING  CYCLE 
THINNING  LEVELS3    INITIAL   -       »»..    *U*SFOUENT  -  80. 
STAND   BEFORE   AND  AFTER  THINNING 


PERIODIC    INTERMEDIATE  CUTS 


BASAL 
AREA 
S4.FT. 

1  19 

57 


107 
78 


97 
115 


AVERAGE 
D. B.H. 
IN. 


6.0 
7.3 


8.3 
9.  1 


11.1 
11.9 


AVERAGE 
HF I8HT 
FT. 

25 

27 


59 
60 


TOTAL 
VOLUME 
CU.FT. 

1190 
630 


1960 
1460 


2140 

2030 


MERCHANT- 
ABLE VOLUME 
CU.FT. 

300 
300 


1680 
1320 


SAUTIHMR 
VetUME 
80. FT. 


2*60 
1910 


iaoo 

1800 


9200 
7500 


TREES 
NO. 


BASAL 
AREA 
SO. FT. 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWTIMBEft 

VOLUME 
BD.FT. 


90 
90 


110 
110 


120 
130 


104 
67 


67 
21 


21 
21 


92 
104 


14.0 
14.8 


16.9 
18.6 


20.2 
21.7 


70 
71 


80 
84 
86 


3400 
2480 


3630 
1410 


1740 
2080 


3230 
2360 


3470 
1360 


1680 
2010 


14400 
10900 


1*000 

7400 


9700 
12200 


TOTAL  YIELDS 

MINIMUM  CUTS  FOR    INCLUSION   IN  TOTAL   YIELDS —     320.   CUBIC   FEET  AND      1500.    BOARD  FEET 
MERCH.   CU.   FT.   -   TREES  6.0   INCHES  D.B.H.    AND  LARGER   TO  4-INCH  TOP • 
BO.   FT.   -   TREES   10.0   INCHES   D.B.H.    AND  LARGER   TO   8-INCH   TOP ■ 


STAND 
AGE 
(YEARS) 

30 
30 


50 
50 


TREES 
NO. 


950 
288 


284 
284 


YIELDS   PER   ACRE   OF   MANAGFD,    EVEN-AGED   STANDS  OF   PONDEROSA  PINE 
SITE    INOEX      70,      20-YEAR  CUTTING  CYCLE 
THINNING  LEVELS=    INITIAL   -       80.,    SUBSEQUENT  -  120. 
STAND   BEFORE   AND  AFTER  THINNING 


PERIODIC    INTERMEDIATE  CUTS 


BASAL 
AREA 
SQ.FT. 

119 

57 


107 
107 


AVFHAGE 
D.B.H. 
IN. 

4>B 
6.0 

7.3 

8.3 


AVERAGE  TOTAL  MERCHANT- 
HEIGHT     VOLUME     ABLE  VOLUME 


25 
27 


45 
45 


1190 
630 


1960 
1960 


300 
300 


1680 
1680 


SAWTIMBER 
VOLUME 
BD.FT. 


1400 
1400 


TREES 
NO. 


BASAL 
AREA 
SO. FT. 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


S  AWT  I MBER 
VOLUME 
BD.FT. 


70 
70 


90 
90 


110 
110 


120 
130 


276 
200 


200 
136 


136 
136 


148 
120 


157 
120 


150 
61 


71 
81 


9.9 
10.5 


12.0 
12.7 


14.2 
15.9 


17.2 
18.4 


58 
58 


64 
69 


B2 
85 


3510 
2930 


4670 
3600 


5070 
2120 


2590 
3070 


3240 
2740 


4390 
3400 


4820 
2020 


2480 
2950 


8400 
8000 


17100 
13800 


22100 
9800 


12900 
16200 


TOTAL  YIELDS 


MINIMUM  CUTS   FOR    INCLUSION   IN  TOTAL   YIELDS —     320.   CUBIC   FEET   AND      1500.   BOARD  FEET 
MERCH.   CU.   FT.   -  TREES  6.0   INCHES   D.B.H.   AND  LARGER   TO  4-INCH  TOP . 
BD.   FT.   -   TREES   10.0   INCHES   D.B.H.    AND  LARGER   TO  E-INCH  TOP. 
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YIELDS   PER   ACRE  OF   MANAGED*    EVEN-AGED  STANDS  OF   PONOEROSA  PINE 


STAND 
AGE 
(YEARS) 

30 
30 


50 
50 


70 
70 


TREES 
NO. 


950 
417 


406 
239 


235 
154 


SITE    INDEX      70,      20-YEAR  CUTTING  CYCLE 
THINNING  LEVELS=    INITIAL   -      110.,    SUBSEQUENT  - 
STAND   BEFORE   AND  AFTER  THINNING 


BASAL 
AREA 
SQ.FT. 

119 
74 

104 

131 


133 
100 


AVERAGE 
D.B.H. 
IN. 

4.8 
5.7 


7.7 

8.5 


10.2 
10.9 


AVERAGE  TOTAL 
HEIGHT  VOLUME 


25 
26 


45 
51 


58 
59 


1190 

800 


2370 
1740 


3220 
2460 


MERCHANT- 
ABLE VOLUME 
CU.FT. 

310 
310 


1900 
1520 


SAWT IMBER 
VOLUME 
BD.FT. 


1400 
1400 


TREES 
NO. 


PERIOOIC    INTERMEDIATE  CUTS 


BASAL 
AREA 
SQ.FT. 


TOTAL  MERCHANT- 
VOLUME     ABLE  VOLUME 
CU.FT.  CU.FT. 


SAWT IMBER 
VOLUME 
BD.FT. 


2990 
2300 


8200 
7300 


90 
90 


110 
110 


120 
130 


154 
104 


104 

32 


32 
32 


133 
100 

114 

129 
49 

58 

68 


12.6 
13.3 


15.1 
16.8 


18.3 
19.  7 


70 
74 


83 
86 


4000 
3050 


4430 
1730 


2150 
2580 


3780 
2890 


4220 
1660 


2060 
2490 


15300 
12200 


20200 
8400 


11200 
14200 


TOTAL  YIELDS 


MINIMUM  CUTS  FOR    INCLUSION   IN   TOTAL  YIELDS —     320.   CUBIC   FEET   AND     1500.    BOARD  FEET 
MERCH.   CU.    FT.   -   TREES  6.0    INCHES   D.B.H.    AND  LARGER   TO  4-INCH  TOP. 
BD.   FT.   -   TREES   10.0   INCHES   D.B.H.   AND  LARGER   TO   8-INCH  TOP. 


STANO 
AGE 
(YEARS) 

30 

30 


50 
50 


TREES 
NO. 


950 
288 


284 
172 


YIELDS   PER   ACRE   OF   MANAGED,    EVEN-AGED   STANDS   OF   PONOEROSA  PINE 
SITE    INDEX      70,      20-YEAR   CUTTING  CYCLE 
THINNING  LEVELS=    INITIAL   -        80.,    SUBSEQUENT  - 
STAND   BEFORE   AND  AFTER  THINNING 


ao. 


PERIODIC    INTERMEDIATE  CUTS 


70 
70 


171 
104 


BASAL 
AREA 
SQ.FT. 

119 

57 


107 
78 


115 

80 


AVERAGE 
D.B.H. 
IN. 


6.0 
7.3 


9.  1 
10.2 


11.1 
11.9 


AVERAGE 
HEIGHT 
FT. 

25 
27 

36 

45 


TOTAL 
VOLUME 
CU.FT. 

1190 
630 


1960 
1460 


2840 
2030 


MERCHANT- 
ABLE VOLUME 
CU.FT  . 

'300 
300 


1680 
1320 


2660 
1910 


SAWT IMBER 
VOLUME 
BD.FT. 


1800 
1800 


9200 
7500 


TREES 
NO. 


BASAL 
AREA 
SQ.FT. 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWT IMBER 
VOLUME 
BD.FT. 


90 
90 


110 
110 


120 
130 


104 
41 


41 
12 


111 

53 

63 

73 

26 

31 
36 


14.0 
15.4 


18.1 
19.8 


21.7 
23.5 


81 
87 


3400 
1670 


2580 
930 


1160 
14  10 


3230 
1590 


2470 
890 


1120 
1360 


14400 
7500 


13400 
5000 


6800 
8700 


TOTAL  YIELDS 


MINIMUM  CUTS  FOR   INCLUSION   IN   TOTAL   YIELDS —      320.   CUBIC   FEET   AND      1500.   BOARD  FEET 
MERCH.   CU.   FT.    -   TREES  6.0    INCHES   D.B.H.    AND  LARGER   TO  4-INCH  TOP. 
BD.    FT.   -   TREES   10.0   INCHES   D.B.H.    AND  LARGER   TO   8-INCH  TOP. 
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YIELDS   PER   ACRE   OF   MANAGED*    EVEN-AGED  STANDS  OF   PONDEROSA  PINE 


STANO 
AGE 
(YEARS) 

30 
30 


50 
50 


TREES 
NO. 


950 
288 


284 
284 


SITE    INDEX     70,      20-YEAR  CUTTING  CYCLE 
THINNING  LEVELS=    INITIAL   -       80.,    SUBSEQUENT  - 
STAND   BEFORE   AND  AFTER  THINNING 


BASAL 
AREA 
SQ.FT. 

119 

57 


107 
107 


AVERAGE  AVERAGE  TOTAL  MERCHANT- 
D.B.H.        HEIGHT     VOLUME      ABLE  VOLUME 


4.8 
6.0 


8.3 
8.3 


25 
27 


36 
45 


1190 
630 


1960 
1960 


300 
300 


1680 
1680 


SAWTIMBER 
VOLUME 
BD.FT. 


1400 
1400 


TREES 
NO. 


PERIODIC    INTERMEDIATE  CUTS 


BASAL 
AREA 
SQ.FT. 


TOTAL  MERCHANT- 
VOLUME     ABLE  VOLUME 


SAWTIMBER 
VOLUME 
BD.FT. 


70 
70 


90 
90 


276 
200 


200 
82 


148 
120 


157 
80 


9.9 
10.5 


12.0 
13.4 


58 
58 


64 
69 


3510 
2930 


4670 
2460 


3240 
2740 


4390 
2330 


8400 
8000 


17100 

9800 


110 
110 


120 
130 


82 
24 


24 
24 


15.5 
17.3 


18.9 
20.4 


78 
80 


3700 
1390 


1730 
2090 


3530 
1330 


1660 
2010 


17100 
6900 


9200 
11700 


TOTAL  YIELDS 


MINIMUM  CUTS   FOR    INCLUSION   IN  TOTAL   YIELDS —     320.   CUBIC   FEET   AND      1500.    BOARO  FEET 
MERCH.   CU.   FT.   -   TREES  6.0   INCHES  D.B.H.   AND  LARGER   TO  4-INCH  TOP. 
BD.   FT.   -   TREES   10. 0   INCHES   O.B.H.    AND  LARGER  TO   8-INCH  TOP . 


STAND 
AGE 
(YEARS) 

30 
30 


TREES 
NO. 


950 
417 


YIELDS   PER   ACRE  OF   MANAGED,    EVEN-AGED   STANDS  OF   PONDEROSA  PINE 
SITE   INDEX      70,      20-YEAR  CUTTING  CYCLE 
THINNING  LEVELS=    INITIAL   -      110.,    SUBSEQUENT   -  100. 


ENTIRE   STAND   BEFORE   AND  AFTER  THINNING 


BASAL 
AREA 
SQ.FT. 

119 

74 


AVERAGE  AVERAGE  TOTAL  MERCHANT- 
D.B.H.        HEIGHT     VOLUME     ABLE  VOLUME 


4.8 
5.7 


25 
26 


1190 
800 


310 
310 


SAWTIMBER 
VOLUME 
BD.FT. 

0 
0 


TREES 
NO. 


PERIODIC   INTERMEDIATE  CUTS 


BASAL 
AREA 
SQ.FT. 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWTIMBER 
VOLUME 
BD.FT. 


50 
50 


70 
70 


406 
239 


235 
154 


131 

94 


133 
100 


7.7 
8.5 


L0.2 
10.9 


45 
51 


58 
59 


2370 
1740 


3220 
2460 


1900 
1520 


2990 
2300 


1400 
1400 


8200 
7300 


90 
90 


110 
110 


120 
130 


154 

62 


62 
19 


19 
19 


133 
66 


90 
34 


41 

47 


12.6 
14.0 


16.3 
18.  I 


19.8 
21.4 


71 
75 


79 
81 


84 
86 


4000 
2050 


3120 
1210 


1510 
1830 


3780 
1940 


2980 
1160 


1460 
1770 


15300 
8500 


15000 
6200 


8300 
10600 


TOTAL  YIELDS 

MINIMUM  CUTS  FOR   INCLUSION   IN  TOTAL   YIELDS —     320.   CUBIC  FEET  AND     1500.   BOARD  FEET 
MERCH.   CU.   FT.   -   TREES  6.0   INCHES  D.B.H.   AND  LARGER   TO  4-INCH  TOP. 
BD.   FT.   -   TREES   10.0   INCHES   D.B.H.    AND  LARGER  TO   8-INCH   TOP . 


21 


APPENDIX  3 
Listing  and  Output  of  Program  PONCHK 


C   TO  OBTAIN  DATA  FOR  COMPUTING  FUNCTIONS  USED   IN  PONYLD. 
C 

C    DEFINITIONS   OF  VARIABLES. 
C 

C  ADM ( I )    =   AVERAGE   DBH  OF  RESIDUAL   STAND  I- 

C  AGE  =   AVERAGE   AGE  OF  MAIN  STAND  TREES  ON  PLOT. 

C  AHT(I)    =    AVERAGE   HEIGHT  OF   RESIDUAL    STAND  I. 

C  AREA    =   AREA  OF  PLOT. 

C  BASH)    =   BASAL   AREA  OF  RESIDUAL   STAND  I. 

C  BDFC(I)   =  H  BD.   FT.   REMOVED  TO  GET  RESERVE  STAND   I.   FROM  TREE  VOLS. 

C  BFO(I>    =   M   BD.    FT.    OF    STAND    I    AS    SUM   OF    TREE  VOLUMES. 

C  BHT(I)   =   AVERAGE  HEIGHT  OF  DOM.   AND  CODOM.   TREES  ON  PLOT. 

C  BIG(l)   =  NUMBER  OF  DOM.    AND  CODOM.   TREES  PER  ACRE. 

C  CC  =  CROHN  CLASS  OF    INDIVIDUAL  TREE. 

C  CFMC(I)    =   MERCH.   CU.   FT.   REMOVEO  TO  GET  RESERVE  STAND   I,    TREE  VOLS. 

C  CFO(I>   =  MERCH.   CU.    FT.  OF   STAND   I    AS  SUM  OF   TREE  VOLUMES. 

C  CRNU)    =  CROHN  CLASS  OF  TREE  I. 

C  DBH  =  D.B.H.    OF    INDIVIDUAL  TREE. 

C  OEN(I)   =   NUMBER  OF  TREES    IN  RESIDUAL   STAND  I. 

C  DIAMII)    =  D.B.H.   OF   TREE  I. 

C  FCBF(I)    =  BD.    FT.    PER  TOTAL  CU.  FT. 

C  FCCM(I)    =  MERCH.   CU.   FT.   PER   TOTAL   CU.  FT. 

C  GRPS(I)    =  NUMBER  OF   TREES   IN  GROUP  I. 

C  HGTUI    =   TOTAL  HEIGHT   OF  TREE  I. 

C  HT  =  HEIGHT  OF    INOIVIOUAL  TREE. 

C  IRNDII)    =  ARRAY  OF   PSEUDORANDOM  NUMBERS. 

C  KCINBIII    =  CUBIC   FEET   PER  M  BD.    FT.,    TREES  REMOVED  ONLY. 

C  NPLT  =   PLOT  NUMBER. 

C  NTIM  =  PAST  OR   PRESENT  CONDITION  OF  PLOT. 

C  NHHN  =  CODED  RECORD  OF  WHETHER  PLOT  HAS   BEEN  THINNED  PREVIOUSLY. 

C  PRET(I)    =   PCT.   OF  ORIGINAL  TREES    IN  RESERVE   STAND  I. 

C  RNCNII)    =  CROHN  CLASSES  OF   STAND  TREES »    IN  RANDOM  ORDER. 

C  RNDTR(I)    =   D.B.H.   OF  STANO  TREES,    IN  RANDOM  ORDER. 

C  RNHT I  I )    =   HEIGHTS   OF    STAND  TREES ,    IN  RANDOM  ORDER. 

C  SUBC(I)    =  MERCH.   CU.   FT.   OF   SUBSAHLOG  TREES  REMOVED  TO  GET  STAND 

C  I ,    AS    SUM  OF    INDIVIDUAL    TREE  VOLUMES. 

C  SUBMFIIJ    =   MERCH.  CU.   FT.   SUBSAHLOG  TREES    IN  STAND   I   AS  A    SUM  OF 

C  INDIVIDUAL    TREE  VOLUMES. 

C  TOCCUI    =   TOTAL  CU.   FT.   REMOVED  TO  GET  RESERVE  STANO   I,    TREE  VOLS. 

C  TOTOU)    =  TOTAL  CU.   FT.  OF  STAND   I    AS   SUM  OF   TREE  VOLUMES. 

C 

DIMENSION  PRE  T I  9 ) ,GRPS<9) »T0CC(9)  ,  SUBMF ( 9 ) »  SUBC ( 9 ) ,CFO ( 9 ) ,CFMC(9)  , 
IBDFC (9)  .BFOI9) , TOT Of  9 ) ,OEN( 9 ) , BAS 1 9 ) , AHT t  9  I , ADM  I  9 ) . B I G ( 9 ) ,BHT(9) ,F 
2CCH(9t , FCBF ( 9 ) , KPRET ( 9 ) , K ADM ( 9 ) , K B AS ( 9  I . KBHT ( 9  ) , K AHT 1 9 ) , K TOT 0 ( 9  I ,K 
3FCCM ( 9 ) ,KFCBf (9)  ,K0ENI9) ,DI AM ( 1000) , HGT I  1000 t .CRN (  1000  I  . I RND (  1000) 
4.RNDTR  (1000)  ,  RNHT  (  1000)  .RNCNl  1000  ) ,KC  JN819) 

C 

C    INITIALIZE    BATCH  VARIABLES. 
C 

GRPSIl)  =  10.0 
GRPS ( 2  )  =  6.0 
GRPSI3)  =  5.0 
GRPSI4)  =  4.0 
GRPS I 5 )  =  3.0 
GRPSI6)  =  2.0 
GRPS ( 7 )  =  3.0 
GRPS  1 8 )  =  4.0 
GRPS19)    =  1.0 

C 

C    INITIALIZE    PLOT  VARIABLES. 
C 

1    TERM  =  21.0 
AREA  =  0.0 
00  3  1=1,1000 
CRN  II)   =  0.0 
DIAMII)    =  0.0 
HGT I  I  )    =  0.0 
I RND ( I  1    -  0 
RNDTRI I  1=0.0 
RNCNII)   =  0.0 
RNHT ( I  )    =  0.0 
3  CONTINUE 
DO  6  1=1,9 
AOM ( I )    =  0.0 
AHT 1  I )  = 
BASH)  = 
BDFC (  I  ) 
BFO(I)  = 
BHT(I)  = 
BIG! I )  = 
CFMCI I ) 
CFO( I )   =  0.0 
DEN  t  I  >    =  0.0 
FCBF  1 1  )    =  0.0 
FCCMt 1 >    =  0.0 
K  ADM (  I )    =  0 
KAHTU)    =  0 
K  B  A  S  (  1  I    -  0 
KBHT ( I )    =  0 
KCINBI I )    =  0 
KDENI I )    =  0 
KFCBF I  I )    =  0 
KFCC  M ( I  )    =  0 
KPRET!  I  )    =  0 
KTOTOII)    =  0 
PRE  T ( I  )    =  O.U 
SUBCU)   =  0.0 
SUBMF I  I  1=0.0 
TOCC (I  1    =  0.0 
TOTOU  1    =  0.0 
6  CONTINUE 

READ    15,10)    NPLT, NTIM, NHHN, AREA, AGE 
10  FORMAT    I I3,212,F5.3,F4.0) 


0.0 
0.0 

=  0.0 
0.0 
0.0 
0.0 

=  0.0 


00   17  L=l,1000 

READ    (5,15)    DBH ,HT  »CC 

15  FORMAT    (F3.1  ,2F3.0) 
IFID8H   .EQ.   99.9)   GO  TO  18 

AA  =  L 

MULT  =  ACRE   *   A  A 

I F ( MULT    .LE.    0)    MULT    =  I 

DO   16  I=KA,HULT 

D I  AMI  I >    =  DBH 

HGT ( I )    =  HT 

CRN (  I )    =  CC 

16  CONTINUE 

KA  =   MULT   +  1 

17  CONTINUE 

18  TREES  =  MULT 

C 

C    COMPUTE    VALUES  OF    STAND   PRIOR    TO  THINNING. 
C 

PRETI9)    =  100.0 
DO  34   1=1, MULT 
SOFT  =  0.0 
TREBF   =  0.0 
TREMCF  =  0.0 
TRETOT   =  0.0 

SOFT  =  0.0054542   *  DIAMII)    *  DIAMII) 
D2H  =    DIAMII)    *   OIAM(I)    «   HGT I  I  ) 

C 

C  IF   STATEMENT   AND  STATEMENTS   FOR  TRETOT   ARE  SP EC  I E S- SP F C I F IC 

C 

IFID2H  ,GT.   6000.0)    GO  TO  20 
TRETOT  =  0.002213  *  D2H  ♦  0.0302B6 

GO  TO  22 

20   TRETOT   =    0.002474   «   D2H  -  1.557103 

C 

C  CHANGE   6.0    IF    OTHER    MINIMUM    D.B.H.    WANTED  FOR   MERCH.    CU.  FT 

C 

22   IFIDIAMII)    .LT.    6.0)    GO  TO  30 

C 

C  if   STATEMENT  AND   STATEMENTS  FOR   TREMCF   ARE  SPECIES-SPECIFIC 

C 

IFID2H   .GT.   6700.0)   GO  TO  24 
TREMCF   =  0.002297  *  D2H  -  1.032297 

GO   TO.  26 

24   TREMCF  =  0.002407   *  D2H  -  2.257724 

C 

C  CHANGE    10.0   IF   OTHER   MINIMUM  D.B.H.   WANTED  FOR   BD.  FT. 

C 

26   IFIDIAMII)    .LT.    10.0)   GO  TO  30 

C 

c  ip   STATEMENT   AND  STATEMENTS   FOR  TREBF   ARE  SPECIES-SPECIFIC. 

C 

IFID2H   .GT.    16000.0)   GO  TO  28 

TREBF   =    (0.012331   »  D2H  -  34.167170)    «  0.001 

GO  TO  30 

28  TREBF  =    (0.016318  *  D2H  -  99.212720)    *  0.001 

C 

C  COMPUTE   ACTUAL   PLCT   VALUES   AS   SUMS  OF    INDIVIDUAL   TREE  VALUES. 


30  AHT I  9 )  =  AHT ( 9 ) 

BAS I  9 )  =  BAS ( 9 ) 

BF0I9)  =  8F0I9) 

CF0I9)  =  CF0I9) 


HGT { I ) 
SOFT 
TREBF 
TREMCF 


I F ( NPLT    .EQ.    888)    GO  TO  400 


T0T0I9)    =   T0T019)    *■  TRETOT 
0ENI9)   =  DEN  1 9 )    +  1.0 
1  F ( CRN ( I )    .GT.   2.0)   GO  TO  32 
BHT (9)    =    BHT (9)    ♦    HGT ( I ) 
BIGI9)   =  BIGI9I    +  1.0 

C 

C  CHANGE   10.0   IF  OTHER  MINIMUM  D.B.H.  USED   FOR  BD.  FT. 

C 

32   IFIDIAMII)    . GE .    10.01    GO  TO  34 

SUBMF(9)    =   SU8MF ( 9 )    +  TREMCF 
34  CONTINUE 

A0M(9)    =  8ASI9)    /    (0.0054542  *  DEN  1 9 )  ) 

IF(ADM(9)    .LE.   0.0)   GO  TO  36 

ADM  I  9)   =  SQRT(ADM(9>> 
36   IF(DEN(9)    .LE.   0.0)   GO  TO  38 

AHT ( 9 )    =   AH  T ( 9  J    /    DEN  I  9 ) 
38   IF(BIG(9)    .LE.   0.0)   GO  TO  50 

BHT  <  9 )    =   8HT19)    /  BIGI9) 

C 

C  COMPUTE   VALUE  OF   MOOULO  FOR  GENERATOR   SUITABLE  FOR   PLOT  DENSITY 

C 

50  DO  55  1=1,12 
TEH  =  I 

THOPR  =   2.0  **  TEM 

IFITWOPR   .GE.   TREES)   GO  TO  60 
55  CONTINUE 
60  NTHO  =    THOPR    *  0.5 

C 

C   GENERATE    PSEUDORANDOM  NUMBERS. 

C 

DO  85   1=1, MULT 
60  NDIV  =    (17.0   *  TERM  ♦   3.0)   /  TWOPR 
NTERM  =  TERM 

NTERH  =   117  *  NTERM  ♦  3)   -  NTWO  *  NDIV 
TERM  =  NTERM 

I F ( NTERM   .LE.   0)    GO  TO  80 

I F (NTERM   .GT.   MULT)   GO  TO  80 


C  REARRANGE  TREES  AT  RANDOM. 
C 

DO  90  I=1,MULT 
NBR  =    I  RND ( I  1 
RNDTRINBR)   =   D I AM ( I ) 
RNHT ( NBR )    =  HGT ( I ) 
RNCN1NBR)    =   CRN(  I  ) 
90  CONTINUE 


SMALLEST  FIRST. 


-  22  - 


00    150  M=l,6 

NOGP  '  TREES   /  CRPS(M) 

ANOGP  ■  NOGP 

TEM  =  ANOGP  «  GRPSIM) 

IF(TEM   . EQ.   TREES )   GO  TO  92 

NOGP   =  NOGP  ♦  1 

ANOGP  =  ANOGP   ♦   I .0 

C 

C   SELECT  LARGEST   TREE    IN  EACH  GROUP. 


92  LX  *  1 

1GP  =  GRPS(M) 
MX  =  IGP 

DO  125  IJ=1,N0GP 
00  94  !L=l.MULT 
JK  =  MULT  -   I L  +  I 
00  94  NK^LX.HX 
NflR    =  NK 

I F ( RNOTR 1 NBR  J  .EQ 
94  CONTINUE 


OIAMI JK> )    GO  TO  98 


C  COMPUTE   BASAL   ARE  A  ANO  VOLUME  OF  EACH  SELECTEO  TREE. 

C 

98    SOFT    =  0.0 
TREBF   =  0.0 
TREHCF  =  0.0 
TRETOT  =  0.0 

SOFT  =  0.0054542   •  RNOTR 1 NBR 1    *  RNOTR ( NBR  ) 
02H  =  RNOTR ( NBR )    •   RNOTR I N8R )    *  RNHT(NBR) 

C 

C  IF   STATEMENT  AND   STATEMENTS  FOR  TRETOT   ARE  SPECIES-SPECIFIC. 

C 

1FID2H   .GT.   6000.01   GO  TO  100 
TRETOT  =  0.002213  *  02H  ♦  0.030288 
GO  TO  102 

100  TRETOT   =   0.002474  •  02H  -  1.557103 


-CHANGE   6.0    IF   OTHER    MINIMUM   D.8.H.   USEO  FOR  MERCH.   CU.  FT. 


102    I F  I  RNDTR ( NBR )    .LT.    6.01    GO  TO  110 
 IF   STATEMENT   AND  STATEMENTS   FOR  TREMCF  ARE   SPEC  I E S- SP EC  I F  IC . 


IF1D2H  .GT .   6700.0)   GO  TO  104 
TREHCF  =  0.002297   •  D2H  -  1.032297 
GO  TO  106 

104  TREMCF  =  0.002407  »  02H  -  2.257724 
 CHANGE   10.0   IF  OTHER  MINIMUM  O.B.H.  USED  FOR  BO.  FT. 

106   I F i RNDTR ( NBR )    .LT.   10.0)   GO  TO  110 

 IF   STATEMENT  AND   STATEMENTS  FOR   TREBF   ARE   SPECIES- SPECIFIC. 

IF(02H   . GT.    16000.01   GO  TO  108 

TREBF   =    10.012331  •  02H  -  34.1671701 
GO  TO  110 

108  TREBF  =    10.016318  *  02H  -  99.212720) 

110  AHTIM)    =   AHTIMI    ♦  RNHT(NBR) 

I F ( RNCN ( NBR)    .GT.  2.0)   GO  TO  115 

BHTIM)    =   BHT(M)    ♦  RNHT(NBR) 

BIG(M)   =  BIG(M)   *  1.0 

115   BAS(M)    =   BASIM)    ♦  SQFT 

TOTO(M)    =  TOTO(H)  ♦  TRETOT 

CF01M)   =  CF0(M1   ♦  TREMCF 

BFO(M)    =  BFOfM)   ♦  TREBF 

DEN(M)    =   DEN ( Ml    ♦  1.0 


0.001 
0.001 


C  CHANGE    10.0  IF  OTHER  MINIMUM  D.B.H.  USED  FOR  8D.  FT. 

C 

IF (RNDTR (NBR 1    .GE.    10. 0)    GO  TO  120 

SUBMFIM)    =   SUBMFIMl   ♦  TREMCF 
120  LX  =  MX  ♦  I 
125  MX  =  MX  ♦  IGP 

C 

C   COMPUTE    ACTUAL   PLOT  VALUES. 
C 

AOM(M)    =   BAS(M)    /    (0.0054542   *  OENIM)) 

I F (ADM  I M )    . LE .   0.0)   GO  TO  130 

ADM(M)    =    SORT ( ADM ( Ml ) 
130   1F(0EN(M)    .LE.   0.0)   GO  TO  135 

AHT(M)    =   AHTIM)    /  DENIM) 
135   IF(BIGIM)    .LE.   0.0)   GO  TO  140 

BHT(M)    =   BHT(M)    /  BIG(H) 
140  PRET(M)    =    ( DEN ( Ml   /  DEN(9)J    *  100.0 

C 

C  COMPUTE   VALUES  FOR  PART  OF   STAND  TO  BE  REMOVEO. 
C 

SUBC(M)    =    SUBMFI9)    -  SUBMFIMl 
CFMC(M)   =  CF0(9)  -  CFO(M) 
BDFC(M)    =   BF0(9)    -  8F0(M> 
TOCC(H)    =    T0T0(9)    -  TOTO(M) 


150  CONTINUE 


ONE   AT   A  TIME. 


DO  200  M  =  7,8 

NOGP  =  TREES  /  GRPS(M) 

ANOGP  =  NOGP 

TEM  =  ANOGP  *  GRPS1M) 

I F I TEM   .EQ.   TREES)   GO  TO  152 

NOGP   =  NOGP  ♦  1 

ANOGP  =  ANOGP  ♦  1.0 

C 

C   SELECT   SMALLEST  TREE    IN  EACH  GROUP. 


152  LX  =  1 

IGP  =  GRPS1M) 
MX  -  IGP 

00  175  IJ=1,N0GP 
DO  154  IL=1,MULT 
DO  154  NK=LX, MX 
NBR  =  NK 

I F (RNOTR ( NBR 1  .EO. 
154  CONTINUE 


OIAMI  ID  )    GO  TO  158 


C  COMPUTE  BASAL  AREA  AND  VOLUME  OF  EACH  SELECTED  TREE. 
C 

158   SQFT  =  0.0 
TREBF   =  0.0 
TREMCF  =  0.0 
TRETOT  =  0.0 
SOFT   =  0.0054542 
02H   =   RNDTR ( NBR ) 


 IF   STATEMENT   AND  STATEMENTS   FOR  TRETOT   ARE  SPECIES-SPECIFIC. 

IFI02H   .GT.   6000.0)   GO  TO  160 
TRETOT  =   0.002213  •  02H  ♦  0.030288 
GO  TO  162 

160  TRETOT   =   0.002474  •  D2H  -  1.557103 
 CHANGE  6.0   IF   OTHER  MINIMUM  O.B.H.   USED  FOR  MERCH.  CU.  FT. 

162   I F (RNOTR I NBR )    .LT.   6.0)   GO  TO  170 

 IF   STATEMENT   ANO  STATEMENTS  FOR  TREMCF   ARE  SPECIES -SPECIFIC. 

IF1D2H   .GT.   6700.0)   GO  TO  164 
TREMCF   =   0.002297  »  02H  -  1.032297 
GO  TO  166 

164   TREMCF   =  0.002407   •  02H  -  2.257724 


-CHANGE    10.0    IF   OTHER   MINIMUM  O.B.H.   USED   FOR   BD.  FT. 


166   1 F I RNDTR ( NBR )    .LT.    10. 0)    GO  TO  170 
 IF    STATEMENT    AND    STATEMENTS    FOR   TREBF    ARE   S P EC  I ES- SP E C  I  F  I C . 


IFI02H 
TREBF  = 
GO  TO  170 
168   TREBF  =  (0.016318 
170  AHTIM)    =    AHT ( M 1  + 
I F I RNCN I NBR)  .GT. 
BHTIM)    =  BHTIM)  ♦ 
BIGIM)   =  BIGIH)  ♦ 
172    BASIM)    =   BASIM)  + 
TOTOIM)    =  TOTQIM) 
CFO(H)    =  CFOIM)  + 
BFOtM)    =   BFOIM)  ♦ 
OEN(M)    =    DEN ( M )  ♦ 


*  02H  -  99.212720) 
RNHT (NBR) 

2.0)   GO  TO  172 
RNHT (NBR ) 
1.0 
SOFT 

♦  TRETOT 
TREMCF 
TREBF 
1.0 


-CHANGE    10.0   IF  OTHER  MINIMUM  D.B.H.   USEO  FOR   BD.  FT. 


AHT I M ) 
BHTIM) 
BIGIM) 


I F I RNDTR I NBR )    .GE.    10. 0)   GO  TO  174 
SUBMFIM)    =   SUBMFIM)   ♦  TREMCF 

174  LX  =  MX  +  I 

175  MX   =  MX   +  IGP 
AHT(T)    *   DEN ( 9 )    -  AHT(M) 
BHTI9)   *   BIG<9>   -  BHT(M) 
BIG(9)   -  BIGIM) 

BASIM)  =  BASI9)  -  8ASIM) 
TOTO(M)  =  T0T0(9)  -  TOTOIM) 
CFOIM)  =  CF0I9)  -  CFOIM) 
BFOIM)  =  BF019)  -  BFOIM) 
DEN(H)  =  DEN ( 9)  -  DEN ( M) 
SUBMFIM)    =    SUBMFI9)    -  SUBMFIM) 

C 

C  COMPUTE    ACTUAL   PLOT  VALUES. 
C 

ADM ( H)    =  BASIM)    /    (0.0054542  *  DEN ( H )  ) 

IFIAOMIM)    .LE.   0.0)   GO  TO  180 

ADM  I M )    =   SQRT I AOM( Ml ) 
180   I F I  DEN  1 M )    .LE.   0.0)   GO  TO  185 

AHT(M)    =   AHTIM)    /  DENIM) 
185   IFIBIGIM)    .LE.   0.0)   GO  TO  190 

BHTIM)    =  BHT(M)    /  BIGIM) 
190   PRE  T  I  M )    =    I  DEN  I M)    /    DENI 9 ) )    *  100.0 

C 

C  COMPUTE  VALUES  FOR  PART  OF  STAND  TO  BE  REMOVED. 
C 

SUBCIM)    =    SUBMF(9)    -  SUBMF(M) 
CFHC(M)    =  CF0I9)   -  CFOIM) 
BDFC IM)    =   BF0(9)    -  BFOIM) 
TOCC(H)   =   T0T0I91    -  TOTO(M) 
200  CONTINUE 


C  COMPUTE   VOLUME  CONVERSION  FACTORS   FOR   EACH  PLOT . 
C 

DO  260  1=1.9 

IF(BOFCU)    .EQ.   0.0)   GO  TO  255 

KCINBII)   =    ((CFMC(I)   -  SUBCII1)    /   BDFCII))   *  10000.0 
255   IF(TOTOII)    . EO.   0.0)   GO  TO  260 

FCCMI I )    =  CFOI I )    /  TOTOI I ) 

FC8FII)   =    (6F0(I)   *   1000.0)   /  TOTO(I) 
260  CONTINUE 

C 

C   SET  NUMBER  OF  DECIMAL   PLACES  TO  BE  RETAINED  FOR  PUNCHING. 


100.0  ♦  0.5 

100. 0  ♦  0.5 

100.0  *  0.5 

100.0  ♦  0.5 


KFCBF ( I )   =  FCBF 1 1  I 

KFCCMI I )    =  FCCMI I ) 

KPRETI  I  )    =  PRETI  1  ) 

KTOTOI I >  =  TOTOI I ) 
2^0  CONTINUE 


1000.0  ♦  0.5 
1000.0  ♦  0.5 
1000.0  +  0.5 
100.0  ♦  0.5 


C   WRITE   RESULTS,   ONE   PAGE   PER  PLOT. 

C 

WRITE  (6.300) 

300  FORMAT    I IH  I , / ,44X ,4 OH VOLUME    ANALYSIS   OF    PONDEROSA   PINE  STANOS) 

HRITE    16,302)    NPLT , NTI M, NWHN 
302  FORMAT    ( 1H0 ,9X , 1 1HPL0T  NUMBER , I4.1X, 1H-, 12. 1X.1H-.I2) 

HRITE  (6,304)  PRET 
304  FORMAT    I  1 HO , / ,  1 X  ,  1 9HPCT .   TREES  RETAINED.4X.9F10. 2) 

WRITE  (6,306)  DEN 
306  FORMAT    (1H0.12HN0.   OF   TRE E S ,  1 1 X ,  9F  1 0 .  0  ) 

WRITE  (6,308)  BIG 
308  FORMAT    (1H0.19HNO.    DOM.   AND  CODOM . , 4X , 9F 1 0 . 0 ) 

WRITE  (6,310)  ADM 
310  FORMAT    I 1H0, 14H AVE RAGE   0. B. H. , 9X , 9 F 10 . 2  I 

HRITE  (6,312)  AHT 
312  FORMAT    (  1H0 , 14HAVERAGE  HE IGHT , 9X , 9F 10 . 2 ) 

HRITE  (6,314)  BHT 
314  FORMAT    IIH0.18HAVE.   HT.   0.   AND  C . , 5X , 9F 10 .2 ) 

HRITE  (6.316)  BAS 
316  FORMAT    ( 1H0, 10HBASAL   ARE A . 1 3X , 9F 1 0 . 2 ) 

WRITE  (6,320) 

320  FORMAT    ( IHO, 35HACTUAL   VOLUMES,   FROM  TREE  VOLUMES  -) 

WRITE  16,322)  TOTO 
322  FORMAT    ( IHO , 5X , 1 3HT0T AL  CU.   FT . , 5X , 9F 10.2 ) 

WRITE  (6,324)  CFO 
324  FORMAT    (  IHO, 5X , 14HMERCH.   CU.   FT. ,4X, 9F  10-2 ) 

WRITE  (6,326)  BFO 
326  FORMAT    (1H0,5X,9HM  BD.   FT . , 9X , 9F 10 . 2  I 
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WRITE    (6.3281  SU8MF 
328  FORMAT   < 1HO, 5X , 1 6HCU.    FT.,    SUBSAWLOG, 9F 10.2 » 

WRITE  (6,3301 
330  FORMAT    ( 1H0, 29HC0NVERS ION  OF   TOTAL   CU.    FT.  -I 

WRITE    (6,332)  FCCM 
332  FORMAT    ( 1H0 , 5X , 16HT0  HERCH  CU.   FT . , 2X , 9F 10. 3 ) 

WRITE    (6,334)  FCBF 
334  FORMAT   ( 1H0, 5X , 1  OHIO  60.   FT. , 8X , 9F 10 . 3 ) 

WRITE  16,345) 

345  FORMAT    { 1H0 , 36HV0LUHES  REMOVEO,    FROM  TREE  VOLUMES 
WRITE    (6,322)  TOCC 
WRITE    (6,324)  CFMC 
WRITE    (6,326)  BOFC 
WRITE    (6,328)  SUBC 


C   PUNCH  RESULTS  FOR    INPUT   TO  REGRESSION  ANALYSIS. 


DO   365  1=1,9 

WRITE    (7,360)    NPLT ,NT I M, NWHN, KPRE T ( I  I ,KADM( I ) , K B AS ( I  I ■ KBMT ( I ) ,K AHT 
L(I> iKDENII ) ,KTOTO( I ) .KCINBI I ) ,KBHT)9) 
360  FORMAT  (13,212,17,15,316,15,217,161 

365  CONTINUE 

00  375  1=1,9 

WRITE    (7,3701    NPL T ,NTI M.NWHN, KPRE T I  I  I , K AOM ( 9  I , KDEN ( I  I , K ADH ( I  I 
370  FORMAT  (13,212,17,3151 
375  CONTINUE 

DO   385  1=1,9 

WRITE    (7,380)   NPLT , NT  I M , NWHN ,KAOM (  I ) , K RAS ( I ) .KFCCM < I ) , KFCBF (I) 
380  FORMAT    ( I  3 ,2  I  2, I  5 , 3 1 6  I 
385  CONTINUE 

GO  TO  1 
400  CALL  EXIT 

END 


PLOT  NUMBER  62 


VOLUME   ANALYSIS   OF   PONDEROSA   PINE  STANDS 


PCT.   TREES  RETAINED  10.00  17. 14  20.00  25.71  34.29  50.00  65.71  74.29  100.00 

NO.   OF   TREES  7.  12.  14.  18.  24.  35.  46.  52.  70. 

NO.   DOM.    AND  CODOM.  7.  12.  14.  18.  24.  35.  46.  52.  66. 

AVERAGE   D.B.H.  19.86  19.17  19.06  18.90  18.29  17.85  17.50  17.56  16.67 

AVERAGE  HEIGHT  86.86  85.17  85.07  85.39  83.87  83.46  83.02  83.02  82.10 

AVE.   HT.   D.   AND  C.  86.86  85.17  85.07  85.39  83.87  83.46  83.02  83.02  82.29 

BASAL  AREA  15.06  24.04  27.73  35.06  43.81  60.84  76.85  87.45  106.09 
ACTUAL  VOLUMES,    FROM   TREE   VOLUMES  - 

TOTAL   CU.   FT.  584.07  913.54  1051.99  1333.21  1638.01  2261.05  2835.88  3227.72  3873.63 

MERCH.   CU.   FT.  563.06  879.89  1013.10  1283.74  1575.82  2173.82  2724.91  3101.68  3716.73 

M   BD.   FT.  3.23  4.96  5.69  7.19  8.67  11.80  14.62  16.66  19.49 

CU.  FT.,  SUBSAWLOG  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00 
CONVERSION  OF    TOTAL   CI).    FT.  - 

TO  MERCH  CU.   FT.  .964  .963  .963  .963  .962  .961  .961  .961  .959 

TO  BD.  FT.  5.530  5.427  5.412  5.395  5.293  5.219  5.156  5.163  5.032 
VOLUMES  REMOVED,    FROM   TREE  VOLUMES  - 

TOTAL  CU.   FT.  3289.56  2960.09  2821.64  2540.42  2235.62  1612.58  1037.75  645.91  0.00 

MERCH.   CU.   FT.  3153.67  2836.84  2703.63  2432.99  2140.91  1542.91  991.82  615.05  0.00 

M  BD.   FT.  16.26  14.54  13.80  12.30  10.82  7.69  4.87  2.83  0.00 

CU.   FT.,    SUBSAHLOG  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00 


Agriculture— CSU,  Ft.  Collins 


-  24  - 


£  s 

si 

<&  g 

Oh 

u  x 
a 

CS  O) 

_,  3 
.3  cs 

B« 

ni 

O)  73 

to  £j 
35  e 


a)  to 

tJJO  <u 
cS 

S  M 

CO  (-i 

3  j,- 

co  5 

2  s 

.2  c 

£3 
u 

O  73 

3  «M 

to 

CD  CO 


>  !-c 

+3  ° 

cs  0. 

"3  6 

si  cd 

H  £ 

73  O 

5  £ 
CO 

~3  cs 
"3 

'3  a 

a> 


~  *i  w 

*H   *~3  H  H-> 


.3     >H  _ 

>o  o 

0)  ^ 

co 


a! 


c  o 


a+J  so 


CS  O 

§3 


t/J  OS 


P  J<!  w 

c 

tO  .  O 
«  Oh-- 
.O  _,  CS 
2  CN  CO 


3  g 

nO  +* 

hS  «S 

SS  ni 

o  > 
Xi 

to  CD 

+J  -O 
CS 

.3  3 

•+->  i-rt 


CD  0) 

3.5 

>  ni 

a)  o  2 

-a  g 

2  S  " 

-  —  - 


Oh 


0) 

s 

0) 

cxo 

cs  _i 

3  * 

2  P 

3  CD 

a;  3 

C  Oh 


G 

~  ni 
ni  cd 


CD  cd 


to  w 
,3  to 

XI  D 
ni  > 

si 

a>  s-. 
m-h  ni 


5.c 

s§ 

>,  to  3 
-  o  « 
m  °  to 

Dh3  Oh 


•-H  tO 

'  3 
1 

3  73 

s 

73  ^ 
%  to 
•g.  73 
"  a> 

M 
■a  cb 

3  0 

CS  cS 

w  £ 


o 

CD 


.2  N  S 

^-  !  S 

a  pi  oj 

CO  Oh 

m  ^  ^ 

^  S  w 

T3  Oh 

a>  ni  a, 

m  o 


to  M 


CS  i-H 

CM 
lO 

a,  to  o 

*  J  £  * 

■o  ?  o  o 

s  « 

2".  CO 
ni  "o 
O 


Onfe 


— 

Si 


o  <J 
cS 


3 


CD  S 
rj  O 


CU  Mh| 

-7-1  CS 


CD 


T3 


X!  O 


°t1  &H  _ 


CS 

cs  -*  *; 


CO  OS 


CD 

"O  ni 

S  S 

rt  cs 

cc  6 

X 

■a 
id 
o 

& 

>> 


al 

Oh 
!h  X 

CD  h 

Oh 

CS  <u 

P  M 
cs 

K« 

CD  T3 
to  o 
^  cs  -i 

m    tO  O 

•a  £  o 

>  o  o 
b  fa  -2 


CD  to 
tXO  CD 
CS 

cs  Q 

s  ■= 

to  2 

2  s 
.2 

n  o 


>  .a  cs 


>>  >. 

0>  cS 
>  Sh 

+3  ° 

«  Oh 


T3  o 


p  CO 


c  o 


•8  ^ 

S-2 


t«  C35 


'  +H  cS  O 

1/3  ti  o 

CD  0 

Sh  -3  - 

to  .  o 

«  Oh-^ 

jO  CS 

CS  ^  -— 

-3  CN  CO 


O 


CO 

CD  to 

51 

2js 

O  jO 
jH  cs 

o  > 
xi 

to  CD 

+J  hO 
CS 


•S 

aj  P  <d 

3  -a  "S 

CS  CD  CD 

o  to 

CS  .O  CD 

«  Ji  s 

CD  C  O 

^2  CS  Sh 


-2 

3  . 

CS  _> 

si 

CD  S- 

H_(  CS 


ca 

CD 

>.  -o 
«  5 


fl  3 
2  ^ 


CD  Oh 

s « 
"3 

Oh 

s 

o 

CD 


n  10 

CD  — 

•°  2 
6  a 

.3  to 
-  c 

CS  T-j 

c 

73  5 

a>  to 
73 

^  CD 
73  cS 

fl  fl 

CS  CS 

J  s 


CO 


h-5  73 
>>  to  0 
«  o  s 

CS   "  to 

s  is 

OhS  Oh 


o 

CD 


T3 

Si 

73  S  C 
fl  Ph  CD 

cs  Oh 

m  5-1  ^ 
«  cd  ca 

73  O, 

a>  cs  J, 

MCI, 

3  -3  cs 

O  cS 
Sh   0)  73 

o  to  0 

«  C3  T 
«  Oh  C 


CD  to 

tXO  CD 
CS 

3  w> 


CD 

!-.  CD 

o  .a 

73  > 
CD  u 
CD  CD 

P  co 

Sh 
Oh 

H_) 
Oh&H 

o  < 
°  Q 

0  P 

I  CS 


to  cs 

2  6 

.2  o 

CD 

«4-H  CD 

O  .73 

a  h 

o  =— 

to 

CD  to 


>>  >. 
2  =3 

+3  O 

Oh 

"3  6 

in  CD 

H-> 

+-» 

CS  g 

73  O 

C3  (h 

Th  ni 
■3  73 

'3  & 


a 


11 


s  § 

£<s 

O  X2 

Sh  CS 
Oh  0 

^  "cS 
o  > 
^3 

to  CD 
*h  ^3 
CO 

43  0 

•—  CS 


•O  73  C 

CS   CD  CD 

^3  C  to 

CS  -3  CD 

cs  -3  O, 


to 


^  w 

73  S 

2  2  « 

cd  3  O 


3 

CD 

2 

CD 

txo 

3 
CS 

2 

t-"  CD 

lH  73 

s  § 

s  a 

■3  to 
3 
.  3 
to  r-1 

CD  Ph 

iS  73 

3 


CD    Oh        73  iS 


3 

CS 


732  • 

i-c    CD  tO 

O  ■-  <U  Oh 

C^N  LEj  I 
!W    H  ^ 


3      ™  •!= 


-O 

5£ 


CD  CD 

3 


3  rn  y  « 


M3  73 


™  -3 
■w  CM  CO 


- 

tO  C33 


B  3 

^^73 

to  3 
O  ni 


to  fl  .3  .  ^ 


35  3 

CD  Sh 


s 


CO  3 


Oh  l 


CD  1/3 
^  CD 

M 

73  CS 

3  3 

ni  ni 

co  2 


o 

CD 

us 


